D'J: 


.J.r  ^..rx;>j..H    93943 


NPS69-83-003 


NAVAL  POSTGRADUATE  8C 

Monterey,  California 


THESIS 

FILMWISE  CONDENSATION  OF  STEAM  ON 
EXTERNALLY -FINNED  HORIZONTAL  TUBES 

by 

William  M.  Poole 

December  1983 

Th 

esis  Advisor:                 P.  J.  Marto 

Approved  for  public  release;  distribution  unlimited 


Prepared  for: 

National  Science  Foundation 
Division  of  Engineering 
Washington,  DC    20550 


T216783 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 

Commodore  Robert  H.  Shumaker  David  A.  Schrady 

Superintendent  Provost 


This  thesis  prepared  in  conjunction  with  research  supported 
in  part  by  National  Science  Foundation,  Division  of  Engineering, 
Washington,  DC  under  MEA82-03567. 

Reproduction  of  all  or  part  of  this  report  is  authorized. 


Released  as  a 
Technical  Report  by: 


Unclassified 


SECURITY  CLASSIFICATION  Of  THIS  PAGE  n»>i#n  Dmtm  Enfrtd) 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 

NPS69-83-003 


2.  GOVT  ACCESSION  NO 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.     RECIPIENT'S  CATALOG  NUMBER 


4.     T\TLE  (and  Subtlllm) 

Filmwise  Condensation  of  Steam  on 
Externally-Finned  Horizontal  Tubes 


S.     TYPE  OF   REPORT   A   PERIOD  COVERED 

Master's  Thesis; 
December  1983 


6.  PERFORMING  ORG.  REPORT  MUMBER 


7.  AUTHORr*; 

William  M.  Poole 


8.  CONTRACT  OR  GRANT  NUMBER^*; 


MEA82-03567 


9.  PERFORMINO  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California   93943 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 


n.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California   93943 

14.    MONITORING  AGENCY  NAME  &   AODRESSC//  dllUrmi  from  CentrolUng  OHIca) 


12.     REPORT   DATE 

December  1983 


13.  NUMBER  OF  PAGES 

85 


15.     SECURITY  CLASS,   (of  thit  report) 


Unclassified 


15«.     DECLASSIFICATION.  DOWNGRADING 
SCHEDULE 


16.     DISTRiauTION  STATEMENT  (ol  thl»  Report) 


Approved  for  public  release;  distribution  unlimited 


17.    DISTRIBUTION  STATEMENT  (ol  (fi«  mbattmct  ittmrud  In  Block  30,  If  dUltronI  /rooi  Rtpori) 


It.     SURPLEMENTARY  NOTES 


19.    KEY  WORDS  (Continuo  on  rovorao  aido  It  noeaaamrr  and  Idanllty  by  block  nianbar) 

Filmwise  Condensation 
Externally-Finned  Condenser  Tubes 
Vacuum  Condensing  Apparatus 
Heat  Transfer 
Instrumented  Tube 


20.     ABSTRACT  (CanHnua  on  rvraraa  aid*  U  nmcoaaary  and  Idanllty  by  block  numbar) 

The  film-condensation  characteristics  of  a  smooth  tube 
and  six  externally- finned  tubes  having  fins  1  mm  high  and 

1  mm  thick,  and  pitches  of  1.5,  2.0,  2.5,  3.0,  5.0,  and 
10.0  mm,  were  experimentally  tested. 

A  smooth  copper  tube  with  an  active  length  of  133.5  mm, 

an  outside  diameter  of  19.05  mm,  and  an  inside  diameter  of 


DO  .  :i:*\,  1473 


EDITION  OF  1  NOV  6S  IS  OBSOLETE 
S/N  0102- LF- 014-6601 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Dala  Bnlarmd' 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wh^n  Dmtm  Enffd) 


12.7  nun  was  first  tested  to  correlate  the  inside  heat- 
transfer  coefficient  using  the  Sieder-Tate  equation.   The 
leading  coefficient  for  this  equation  was  found  to  be 
0.034  ±  0.001,  and  was  used  to  derive  the  external  con- 
densing coefficient  for  all  of  the  tubes  by  subtracting 
the  inside  and  wall  resistances  from  the  measured  overall 
resistance.   The  condensing  coefficient  was  measured, 
both  at  atmospheric  pressure  and  vacuum  (84  mm  Hg) ,  with 
the  heat  flux  as  a  variable. 

Condensation  data  taken  for  the  smooth  tube  were  com- 
pared with  data  in  the  literature  to  check  the  reliability 
of  the  apparatus  and  the  data-reduction  procedures.   The 
data  for  the  finned  tubes  showed  an  optimum  pitch  of  2.5  mm 


S'N  0102-  LF-  014-  6601 

2  Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PA0C(1Wi*n  Dmtm  Bnfrmd) 


Approved  for   public   release;    disrribu-ion   unlimited 


Filmtfise  Condensation  of   Steam   on 
Externally-finned   Horizontal   Tubes 


by 


William  M.  Poola 

Lieutenant,  United  Scatas  Navy 

B.S.,  Marine  Engi.iearing 

(J.S.  Naval  Academy,  1978 


Submitted  in  partial  fulfillment:  of  ths 
requirements  for  the  deqrae  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 


from  the 


NAVAL    POSTGRADUATE    SCHOOL 
December    1983 


csfS 


\  ■ 


M3 


ABSTRACT 

The  film-aondsnsatioa  characterisuics  of  a  smooth  xube 
and  six  extsrnally-f inned  tubes  having  fins  1  mn  high  and  1 
mm  thick,  and  pitches  of  1.5,  2.0,  2.5,  3.0,  5.3,  and  10.0 
mm,    were   experimentally  tested. 

A  smccth  copper  tube  with  an  active  length  of  133.5  ffim, 
an  outside  diameter  of  19.05  mm,  and  an  inside  diameter  of 
12.7  mm  was  first  tested  to  correlate  the  iiside  hc'at- 
transfer  coefficient  using  the  Sisder-Tare  equation.  The 
leading  coefficient  for  this  equatiDn  was  found  to  be  0.034 
±  0.001,  and  was  used  to  derive  tha  external  condensing 
coefficient  for  all  of  the  tubes  by  subtracting  the  inside 
and  wall  resistances  from  the  measured  overall  resistance. 
The  condensing  coefficient  was  measjrad,  both  at  atmospheric 
pressure  and  vacuum  (84  mm  Hg)  ,  with  the  heat  flux  as  a 
variable. 

Condensation  data  taken  for  the  smooth  tube  were 
compared  with  data  in  the  literature  to  check  the  reli- 
ability of  the  apparatus  and  the  data-reduction  procedures. 
The  data  for  the  finned  tubes  showsi  an  optimum  pitch  of  2.5 
mm. 
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I.    INTRODOCTIDN 

A.       BACKGROUND 

As  combat  systems  aboard  naval  surface  ships  continue  to 
grow  in  importance,  size,  and  weight,  the  feasibility  of 
their  installation  aboard  lightweight  vessels  can  only  be 
realized  by  reducing  other  major  shipboard  weight  require- 
ments. A  significant  achievement  in  reducing  main  propulsion 
weight  was  reached  with  the  advent  of  gas-turbine  propulsion 
plants,  of  which  the  majority  of  our  new  combatants  will  be 
powered.  However,  for  the  remaining  capital  ships  and  subma- 
rines to  be  constructed  using  either  conventional  or  nuclear 
steam  systems,  power  plant  weight  reduction  can  only  be 
realized  by  increasing  the  effectiveness  of  the  individual 
plant  components  and  thereby  reducing  their  size.  Plans  call 
for  future  installation  of  Ranlcine-rycle  waste-heat  recovery 
systems  even  on  those  ships  with  gas-turbine  propulsion,  so 
significant  weight  reduction  could  also  be  realized  there 
through   the   use    of    more   effective   systems. 

Due  to  the  nature  Df  the  condensing  process,  the 
greatest  thermal  resistance  to  condenser-tube  heat  transfer 
occurs  on  the  tube  side,  and  a  thorough  review  of  tube-side 
enhancement  is  presented  by  Berglss  and  Jensen  [Ref-  1  ]- 
Internal  enhancement  of  the  tube  without  regard  to  the 
shell-side  resistance  problem,  however,  would  be  a  wasted 
effort  if  improvements  ware  made  to  the  point  where  the 
external  (shell-side)  resistance  became  the  controlling 
factor  in  the  transfer  of  heat.  Recent  investigations  into 
the  shell-side  condensing  process  and  techniques  to  improve 
this  process  are  thoroughly  reviewed  by  Marto  [Ref-  2]. 
Outside      enhaacement     techniques      include     the     use      of      low 
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integral  fins,  roped  tabes,  fluted  tubes,  and  appli^sd 
coatings   to   promore    drcpwise  condsnsarion. 

Ongoing  research  is  underway  at  the  Naval  Postgraduate 
School  to  further  study  the  shell-side  enhancement  of  marine 
condensers.  An   endurance      apparatus      is      in  operation     to 

examine  the  drop  wise-promoT.ing  effectiveness  of  various 
metallic  and  polymer  coatings  over  prolonged  periods  of 
time.  Reilly  [  Hef -  3]  conducted  a  study  of  imprcvemsnts 
using  various  spirally  fluted  tubes  under  filmwise  condensa- 
tion, and  the  effects  of  tube  bunile  inundation  have  also 
been    studied  by    Kanakis   [Ref.    4]. 

A  test  apparatus  has  been  constructed  by  Krohn  [Ref.  5] 
to  systematically  study  steam  condensation  on  a  single  hori- 
zontal tube.  The  apparatus  is  modelad  after  a  similar  design 
constructed  by  Rose  [Bef.  6],  but  can  also  operate  under  the 
high-vacuum  conditions  found  in  marine  condensers  and 
without    the      presence  of      noncondensable    gases.  This   test 

apparatus  was  instrumented  and  tested  by  Graber  [Ref-  7], 
who  determined  the  Sieder-Tate  coefficient  for  the  test  tube 
length   to   be  0. 029. 

In  experiments  using  low  integral  fins  at  Queen  Mary 
College  of  London,  Rose,  et.  al.  ,  [Ref,  6]  obtained  an 
optimum  pitch  of  2.0  mm  (0-08  in.)  for  a  gi.ven  fin  size  of  1 
mm  (0.04  in)  high  and  0.5  mm  (0.02  in)  wide  for  tests 
conducted  at  one  atmosphere.  With  this  pitch,  he  reported 
vapor-side  heat-transfer  enhancements  of  400^  aad  300%  for 
heat  fluxes  of  0.3  M W/e2  (9.51x10*  Btu/hr-fts)  and  0.8  MW/m2 
(2.54x105  3tu/hr-ft2),  respectively,  and  for  a  vapor 
velocity  of  0.7  m/s  (2.3  ft/sec).  Further  testing  is  needed 
to  compare  these  results  with  similar  fin  geometries  under 
vacuum  conditions,  while  varying  fin  dimensions  (height  and 
width)    in   addition    to  the   fin    pitch. 
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B.       OBJECTIVE 

The   main  objectives   of   this   resaarch   effort   were,    ther' 
fore,    to: 

1)  ensure  a  vacuum -tight  apparatus  so  that  data  could  be 
taken  at  both  atmospheric  and  vacuum  conditions  with  no 
detrimental  effects  due  to  the  presence  of  nonconden- 
sable    gases; 

2)  take  data  for  an  instrumented  smooth  rube  to  verify 
the   Siader-rate  coefficient   obtained   by   Graber; 

3)  take  data  for  a  smooth  tube  to  check  -he  reliabilty 
of  the  apparatus  and  the  data-reduction  procedures  used; 
and 

U)  take  data  for  six  externally-enhanced  tubes  of 
various  fin  pitches  to  obtain  the  relative  op-imum  pitch 
for   a    fixed   fin   gecme-ry. 
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II.     DESCRIPTION    OF   iPPMAIUS 

A-   SISTEH  OVEBVIEW 

The  apparatus  used  for  this  research  was  essentially  th'? 
same  used  in  references  5  and  7,  with  several  noted  modifi- 
cations. A  schematic  sketch  of  the  system  is  shown  in 
Figure  2.1.  Steam  was  generated  in  a  30U.8-m(n  (12-in)  diam- 
eter Pyrex  glass  boiler  by  ten  UOOO-watt,  U80-volt  Watlow 
immersion  heaters.  Passing  througa  a  304.8-mm  (12-in)  to 
152.U-mm  (6-in)  reducing  section,  the  steam  travelled  upward 
through  a  Pyrex  section  2.4U  m  (8.D  ft)  in  length,  around  a 
180-degree  bead,  and  back  down  a  straightening  section  1.52 
m  (5.0  ft)  in  length  before  entering  the  stainless-steel 
test  section.  The  condenser  tube  to  be  tested  was  mounted 
horizontally  in  the  test  section  banind  a  viewport  to  permit 
visual   observation    of  the   condensing    process. 

Steam  that  did  not  condense  on  the  test  tube  passed  into 
a  stainless  steel  auxiliary  condenser,  and  all  condensate 
was  returned  via  gravity  to  the  boiler.  The  auxiliary 
condenser  was  constructed  of  two  9.5-mm  (3/8-in)  water- 
cooled  copper  lines  helically  coiled  to  a  height  of  457  mm 
(18    in)  . 

Cooling  water  for  the  test  tube  was  provided  by  two 
centrifugal  pumps  connected  in  series.  The  water  could  be 
throttled  from  zero  flow  to  0.69  1/s  (11  gpra).  The  maximum 
water  velocity  which  could  be  obtained  through  the  tube  was 
5.48  m/s  (18  ft/sec)  .  A  continuous  supply  of  tap  water  was 
used  for  cooling  the  auxiliary  condenser.  Throttling  the 
flow  of  tap  water  through  the  condsnser  was  the  means  used 
to  vary  the  internal  pressure  of  the  test  apparatus.  The 
water      flow   through      both   the      test   tube      and  the      auxiliary 
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condenser  was  regalated  by  19,1  mm  (3/4  in)  dianeter  needle 
valves  and  measured  by  rotameters  with  full-scale  ranges  of 
0-69    1/s    (1 1    gpm). 

An  air  ejector  provided  for  rBmoval  of  noa condensable 
gases  from  the  auxiliary  condenser  through  a  12.7-mm 
(1/2-in)  line.  The  source  for  ths  air  ejector  was  1.1-MPa 
(160-psig)    house-air    supply. 

B.       SYSTEM    INSTHOMENTATION 

The  input  voltage  through  the  hsatars  was  varied  through 
a  panel-mounted  potentiometer.  a'40-VAC  line  voltage  was 
reduced  by  a  factor  of  100  when  fed  into  a  differen-ial 
input  precession  voltage  attenuator.  The  stepped-dcwn 
voltage  passed  through  a  True-Roo t-Jlean-Sguare  ccnver-er 
stage  on  which  the  integrated  period  was  reduced  to  abou-l:  1 
ms.  The  output  of  the  TRMS  converter  was  then  buffed  and 
compared  to  a  reference  voltage  from  the  potentiometer.  The 
comparator  output  was  fed  to  the  control  input  of  a  Halmar 
silicon-controlled  rectifier  power  supply  which  applied  the 
actual  voltage  to  the  heaters.  The  TRMS  converter  output  was 
also  paralleled  to  a  filter  and  then  input  -o  the  data 
acquisi-ion  system.  This  input,  was  proportional  to  the  p'ower 
supply  output.  A  diagram  of  the  sys-em  is  shown  m  Figure 
2.2. 

The  internal  pressure  of  the  system  was  measured  manu- 
ally by  a  U-tube,  mercury -in-glass  manometer  graduated  in 
millimeters.  Unavoidably,  steam  could  condense  in  the  manom- 
e-csr.  Therefore,  the  varying  height  of  the  water  column  in 
the  mancmeter  needed  to  be  accouniiad  for  when  measuring  the 
system   pressure. 

Temperatures  throughout  the  system  were  leasured  by 
copper-constantan  thermocouples:  six  for  the  wall  of  a 
specially-constructed  test    tube,    two    for   the   steam,      and   one 
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each  for  the  cooling-watar  inlet,  condensate  return,  and 
ambient.  The  calibration  procedure  for  these  thermocouples 
is  described  in  Appendix  A.  The  temperature  rise  through  th9 
tes-  tube  was  measured  by  a  Hewlett-Packard  (HP)  280UA 
quartz   thermometer. 

All  temperature  measurements  were  fed  directly  into  the 
data-acquisition   system   as    described    below. 

C.  DATA    ACQUISITION 

An  HP  3U97A  Data  Acquisition/Control  Unit  was  used  to 
monitor  system  temperaturss .  This  was  interfaced  with  an  HP 
9826A  computer  which  served  as  a  controlling  unit  through  an 
interactive  data-reduction  program  and  user  keyboard 
prompts.  Raw  data  gathered  by  the  data- acquisition  system 
were  stored  on  computer  disks  for  later  reduction  and  evalu- 
ation. 

D.  SYSTEH    HODIFICATIONS 

1 .      Boiler 

The  fiberglass  insulation  was  removed  from  the 
boiler  to  allow  the  operator  to  aore  easily  monitor  the 
water  level.  Although  a  closed-system  design  was  used,  it 
was  still  possible  for  steam  to  escape  via  the  air  ejector 
or  through  the  relief  valve  [Fig.  2.1].  Calculations  showed 
that  the  additional  heat  loss  due  to  the  removal  of  the 
insulation  was  minimal  [App.  3],  and  the  author  felt  this 
loss  was  much  more  acceptable  than  risking  damage  to  the 
immersion  heaters  through  a  low-water  casualty  in  the 
boiler. 
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2-      Con  den  sat  a   ?iEin£ 

As  originally  designed,  draining  rhe  system  required 
brsaking  down  the  condensate  piping.  While  this  was  not  a 
daily  occurrence,  the  procedurs  was  inconvenient  and  there 
existed  also  the  possibility  of  losing  the  vacuum  integrity 
of   the   system   each    tiie    ir    was  done. 

To  avoid  these  problems,  the  existing  fill-line 
valves      were     rearranged      as   shown     in      Figure      2.1.  This 

arrangement  also  added  two  additional  features  to  the 
system: 

1)  the  fill/drain  valve  could  be  opened  during  operation 
to  drain  any  heavy  particulate  matter  from  the  system  - 
similar  to  the  "bottom  blow"  procedure  used  on  Naval 
boilers;    and 

2)  after  extended  periods  of  inactivity  while  opened  to 
the  atmosphere,  the  entire  system  could  be  given  a  thor- 
ough steam-cleaning  by  following  the  procedures  outlined 
in    Appendix   C. 

3.  Vent  Valve 

The  modification  of  the  condensate  return  piping 
necessitated  the  addition  of  a  vent  valve  for  use  when 
filling  or  draining  the  system.  A  U,3-mm  (0.17-in)  needle 
valve  was  installed  on  the  101.6-im  (4-in)  flange  of  the 
test  section.  This  valve  would  also  serve  as  the  tap  for  the 
proposed  sampling  of  noncondensabl=  gas  concentrations  in 
the    system. 

4.  Manometer  Line 

The  original     system   design      used   a      6.4-mra    (1/U-in) 
stainless   steel    tube      angled   down   tD    the      mercury    manometer. 
During   this   thesis,    the    manometer   was   raised   to    eye    level  to 
facilitate      easier    and  quicker      reading.      Replacement   of   the 


19 


stainless-steal  line  by  a  12-7-ini  (1/2-in)  copper  tubi 
reduced  the  possibilty  of  error  caused  by  water  slugs 
building  up  in  -he  smaller  diameter  tube.  The  mere  workable 
copper  was  chosen  over  stainless  steal  to  reduce  the  stiff- 
ness of  the  connecting  line.  This  was  necessary  to 
eliminate  leakage  in  this  part  of  the  apparatus  as  explained 
in    Section    II.  E. 

5-      Pressure   Transducer 

As  an  alternative  to  tha  manometer,  a  Celesco 
strain-gage  pressure  transducer  was  installed  on  the  tesr 
section  flange  next  to  the  vent  valva.  The  calibration  line 
for  the  transducer  is  shown  in  Figure  2.3.  The  author  felt, 
however,  that  the  reliabili-^y  of  tdis  measur aman t  would  not 
be  high  enough  until  a  second,  mora  accurate  transducer  was 
installed.  Once  incorporated  into  tha  system,  though,  these 
transducers  would  provide  automatic  input  to  the  data  acqui- 
sition system,  eliminating  the  requirement  to  maiually  enter 
the    manometer   reading  into    the   data-reduction   program. 

6.      Relief    Valve 

As  originally  designed,  a  5895-Pa  gage  pressure 
(1-0-psig)  relief  valva  was  installed  beneath  a  1.0-m 
(39.U-in)  length  of  12.7-mm  (1/2-in)  srainless-st eel  piping. 
Steam  which  condensed  and  became  trapped  in  this  piping 
would  open  the  valve  at  a  water-cclimn  height  of  only  0.70  ra 
(27.6  in).  Once  opened,  a  back-pcassure  of  0.14  MPa  gage 
(20  psig)  was  required  to  reseat  tha  valve  -  something  unob- 
tainable even  with  an  absolute  vacuum  on  the  inlet  side  of 
the  valve.  To  avoid  this  problem,  tha  valve  was  raised  in 
tha  line  to  a  point  only  76  mm  (3.0  in)  below  the  outlet 
from  the  auxiliary  condenser  section  and  was  replaced  by 
another  6895  Pa  gage  pressure  (1.0-psig)  relief  valve  which 
reseated  at   only   O.OU-MPa   gage    (6-psig)     back- pressure. 
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7-      St_f am-t smoeraturj  Probe 

The  steam  temperature  prooe  was  located  directly 
above  the  test  tube  and  shared  ths  same  inlet  re  the  test 
section  as  the  manometer  line.  When  drawing  a  high  vacuum, 
this  arrangement  allowed  water  which  had  collected  in  the 
manometer  to  be  drawn  bade  into  -he  test  section  where  it 
would  flow  down  the  probe  and  onto  the  tube.  To  prevent  this 
water  and  any  contaminants  picked  up  in  the  manometer  from 
being  deposited  on  the  tube,  the  probe  was  bent  so  that  it 
was    offset    from    the    center    of    the    tast   tube. 

8 •      Cooling    Wa t e r   S^st e m 

A  second  centrifugal  pump  was  added  in  series  to  the 
one  already  installed.  This  pump  boDS-ed  the  maximum  cooling 
water  velocity  through  the  test  tube  from  3.96  m/s  (13 
fr/sec)    to    5.U8    m/s     (18    ft/sec). 

The  10-55  kw  (3  Ton)  air-conditioning  anir  used  with 
the  cooling  water  system  would  energize  at  a  water  tempera- 
ture of  17  "C  (62.6  "  F)  and  secure  when  the  -amperature  was 
reduced  to  13.4  "C  (56.1  "F).  Therefore,  for  a  given  steam 
temperatura  of  50  °C  (122  ^F)  around  the  tube,  the  log- 
mean-temperature  difference  would  vary  by  as  much  as  11%. 
But  the  measured  temperature  rise  of  the  water  through  the 
tube  showsd  very  little  change.  To  avoid  this  -ransient 
problem,  the  air-conditioning  unit  was  not  used,  and  instead 
fresh  tap  water  was  continuously  fed  to  the  saip  while  an 
equal  amount  of  water  was  being  drained,  maintaining  a 
constant  sump  level.  This  method  provided  a  constant- 
temperature  supply  of  cooling  water  to  the  inlet  side  of  the 
test    tube. 

It  should  be  pointed  out  that,  while  the  duty  cycle 
of  the  air-conditioning  unit  was  a  function  of  the  ther- 
mostat  used,      a    more    sensitive     thermostat    would    require   the 
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US9  of  a  more  expensive  cooling  dsvica  than  a  commercial  air 
conditioner.  An  alternative  solution  would  be  to  install  a 
much    larger   sump   in    the  system. 

9 .      Thermopile 

Since  the  temperature  rise  of  the  cooling  water 
through  the  tube  was  a  critical  measurement  in  the  experi- 
ment, a  10-junction  thermopile  was  added  to  neasure  this 
temperature  rise  in  addition  to  the  quartz  thermometer.  As 
will  be  explained  in  Section  V.A.3,  however,  problems  arose 
with  the  thermopile,  and  it  couli  not  be  accurately  used 
during  this   thesis. 

E.       VACODH    INTEGRITY 

One  of  the  objectives  of  this  thesis  was  to  ensure  a 
vacuum-tight  test  apparatus  to  eliminate  the  presence  of  any" 
noncondensable  gases  and  their  detrimental  effects  on  the 
condensing    process,  A   standard   of      no    more   than      a   5.0-mm 

mercury  (0.10-psi)  loss  over  a  24-hr  period  was  considered 
to  be  an  acceptable  tolerance,  bat  obtaining  a  leak  rate 
within  this  tolerance  proved  to  be  the  most  time-consuming 
effort   during  the  research. 

Due  to  the  construction  of  the  apparatus  and  nature  of 
the  experiment,  most  leak-detection  methods  could  not  be 
used.  A  sealing  substance  could  not  be  used  without  risking 
contamination  of  the  interior  of  the  apparatus  which  might 
prohibit    filmwise  condensation    on   tie    test    tube. 

Initially  the  system  was  pressurized  and  the  standard 
soap-solution  test  was  used  to  locate  leaks.  Once  pressur- 
ized, a  liquid-soap  solution  was  applied  to  each  external 
joint  or  fitting  where  a  leak  couli  be  present.  The  higher 
pressure  air  inside  the  apparatus  would  escape  through  any 
leaks      and      produce      bubbles      on        the      applied      soap      film. 
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However,  the  maximum  pressure  which  the  Pyrex  glass  members 
could  tolerate  was  orly  0.07U  MPa  gage  (10.7  psig) ,  and  for 
safety  reasons  the  author  chose  not  to  pressurize  the  system 
to  more  than  0.034  MEa  gage  (5.0  psig)  -  a  pressure  differ- 
ence between  the  atmosphere  and  ths  apparatus  of  only  O.OUO 
MPa     (5.8   psi).  The  numerous   sxtacnal   valves      and    fi-tings 

prohibited  the  use  of  an  evacuated  hood  to  achieve  a  greater 
pressure   difference. 

A  similar  test  could  not  be  used  to  locate  any  vacuum 
leaks  which  were  not  present  when  the  system  was  under  a 
positive  pressure,  as  the  apparatus  was  not  large  enough  to 
permit  the  application  and  observation  of  a  soap  solution  on 
the    interior. 

In  another  attempt  to  locate  the  leaks,  a  National 
Research  Corporation  (NRC)  10  1.6-mm  (4-in)  vacuum  pumping 
system  was  connected  to  the  apparatus.  This  pumping  system 
included  a  Welch  model  1376M  mechanical  pump  and  a  model 
NHS-U  diffusion  pump.  An  NRC  model  521  thermocouple  gage  was 
also  connected  to  the  test  apparatus  and  the  entire  appa- 
ratus was  evacuated  to  0,21  torr  (4.1x10"3  psia) .  Acetone 
was  sprayed  around  all  flanges,  fittings,  and  joints.  A  leak 
around  any  of  these  should  have  produced  a  rapid  rise  in  the 
thermocouple  reading,  but  this  method  also  proved  ineffec- 
tive, probably  due  to  the  large  size  of  the  apparatus 
resulting  in  too  great  a  mean  free  path  for  the  acetone 
molecules   to  travel    from  the   leak    to    the   thermocouple. 

The  next  alternative  was  to  break  the  system  apart  into 
three  main  sections:  the  glass  boiler  and  steam  piping,  the 
stainless-steel  test  section  and  auxiliary  condenser,  and 
the  condensate  return  piping.  The  glass  section  was  blank- 
flanged  and  evacuated  to  an  absolute  pressure  of  0.033  torr 
(6.4x10  -♦  psia) .  The  rate-of-rise  measurement  for  this 
section  showed  a  loss  of  only  0.43  mm  of  mercury  (0.01  psi) 
over  a  30-hr  period.  It  was,  therefore,  concluded  that  any 
leak    in    the   assembled  apparatus    was   not    from    this    section. 
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Once  reniDved,  the  test  section  and  the  auxiliary 
condenser  were  blank-flanged,  pressurized  to  0.10  MPa  gage 
(15  psig) ,  and  immersed  in  a  large  plexiglas  tank  filled 
with  water.  This  test  easily  revealed  a  number  of  small 
leaks  abcut  the  inlet  side  of  rhe  tesu  tube  and  also  around 
the  plug  which  connected  the  condensate  return  piping  to  -he 
base  of  the  auxiliary  condenser.  Replacing  an  0-ring  in  the 
tube  fitting  and  silver  brazing  the  plug  into  place  elimi- 
nated  these   leaks. 

The  same  immersion  test  revealed  small  leaks  in  -iihe 
joints  of  the  condensate  return  piping.  These  leaks  were 
eliminated  by  replacing  all  stainless-steel  ferrules  in  the 
Swagelck    fi-^tings   with  teflon    ferrules. 

Once  reassembled,  considerable  leakage  was  still  indi- 
cated by  a  substantial  overnighi:  rise  in  the  manometer 
level.  The  author  felt  confident  that  this  leak  was  not  in 
the  main  assembly  of  the  apparatus,  but  was  instead  in  the 
manometer  assembly  itself  since  it  could  not  be  pressurized 
for  testing.  As  mentioned  in  section  II.C.U,  the  stainless- 
steel  line  leading  to  the  manometer  was  at  this  zime 
replaced  by  a  12.7-mm  (1/2-in)  soft  copper  tube.  This  tube 
eliminated  the  need  for  two  90-degree  elbows  and  three 
lengths  of  stainless  steel  tubing  ia  the  line  -  an  assembly 
which  proved  too  rigid  to  allow  even  the  slightest  misalign- 
ment   into   the   manometer. 

Open  completion  cf  the  installation  of  this  assembly, 
the  system  was  evacuated  to  an  absolute  pressure  of  92.5  mm 
Hg  and  ever  a  24-hr  period  the  mercury  level  rose  to  only 
94.0  mm.  This  leak  rate  was  well  within  the  acceptable 
tolerance. 
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F.       TOBES   TESrED 

1  •      Instrumented    Tube 

An  instrumented  tube  was  fabricated  from  a  thick- 
walled  copper  tube  with  an  inner  diameter  of  12.70  .urn  (1/2 
in)  and  an  outer  diameter  of  19.05  mm  (3/4  in).  The  tube 
was  cut  into  three  sections  in-o  which  six  holes  were 
drilled  axially  along  the  walls  at  equal  spacings  60  apart. 
These  passages  were  fitted  with  0.094-mm  (3/32-in)  OD  capil- 
laries [Fig-  2.U]  which  were  silver-soldered  into  place,  and 
the  three  sections  of  tube  were  then  soldered  back  into  one 
piece.  Thermocouples      were      fittsi        into      the      capillary 

sections   to   measure    an   average   wall   temperature. 

By  knowing  the  average  wall  temperature,  the  Nusselt 
number  for  the  inside  could  be  computed.  By  computing  the 
gradient  of  the  Nusselt  number  against  the  Sieder-Tate 
parameter,  the  inside  coefficient  could  be  obtained  as  the 
inverse  of  the  gradient.  Figure  2.5  shows  a  photograph  of 
the    instrumented   tube   with    the   installed   thermocouples. 
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2.      Smooth    Tuba 

A  smoDth  tube  with  no  wall  thermocouples  was  also 
tested  to  obtain  the  inside  heaTi-transf  er  coefficient 
through  the  use  of  a  modified  Wilson  Plot  [Ref.  8]. 
Determination  of  the  inside  heax-transf er  coefficient  was 
critical  to  the  experiment,  as  it  was  used  to  obtain  i-he 
outside  condensing  coefficient  for  the  smooth  tube  and  all 
of   the   enhanced    tubes. 

3 •      Finned    T ubes 

To  fulfill  the  main  objective  of  this  thesis,  a 
series  of  six  finned  tubes  was  also  tested  [Fig.  2.6]. 
These  tubes  had  the  same  overall  dimensions  as  those  above, 
but  were  enhanced  with  radial  fins  1  mm  (O.Qi^  in|  high  and  1 
mm  (0.04  in)  thick.  Each  tube  had  a  different  fin  pitch  and 
was  tested  xo  determine  a  relative  optimum  pitch.  Fin 
pitches   tested    were    1.5,    2.0,    2.5,    3,0,    5.0,    and    10.0    mm. 

G.       SYSTEH    OPERATION 

The  tube  to  be  tested  was  cleaned  in  a  warm  solution  of 
Sparkleen  and  then  rinsed  with  tap  water,  which  produced  a 
contaminant-free,  wetted  surface.  The  tube  was  then 
installed  in  the  test  section,  care  being  taken  not  to  touch 
or   contaminate    the    condensing    surface. 

The  system  was  brought  to  operating  pressure  by 
following  the  procedures  of  Appendix  D,  and  data  collection 
began        when        steady-state  conditions        were        achieved. 

Steady-state  conditions  were  determined  by  observing  the 
steam  temperature  measured  by  the  respective  thermocouples. 
When  their  output  voltage  on  the  HP  3497A  reached  a  constant 
value  with  fluctuations  of  only  on2  or  two  microvolts,  it 
was  assumed  that  steady-state  conditions  existed  in  the  test 
apparatus. 
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Data  sets  were  taken  starting  with  a  test-tube  cooling 
water  flow  rate  of  90%  (which  corresponded  to  a  warer 
velocity  of  4.95  m/s  or  16.2  ft/sec) ,  ranging  downward  in 
decrements  of  10%  through  a  minimam  flow  of  20S  (1.16  ra/s, 
3.8  ft/sec),  and  then  upward  from  25%  to  85%  in  increments 
of  10%.  After  adjusting  the  flow  rate,  the  temperature  rise 
through  the  tube  was  monitored  by  observing  the  digital 
output    of      the    quartz  thermometer.  ^/Jhen    this      rise    became 

constant,  another  set  of  data  could  be  taken.  During  data 
runs,  a  slight  rise  in  pressure  accompanied  the  decrease  in 
the  cooling  water  flow  through  nhe  tube,  a  result  of  the 
reduced  heat  flux.  Similarly,  increasing  the  flow  through 
the  tube  caused  a  slight  decrease  in  the  system  pressure. 
This  variance  could  be  anticipated,  and  since  data  were  to 
be  collected  at  a  constant  pressure,  it  was  easily  compen- 
sated for  by  throttling  the  flew  of  cooling  water  through 
the  auxiliary  condenser  one  or  two  percentage  points  on  the 
rotameter. 

Something  which  could  not  be  anticipated,  however,  were 
sudden  fluctuations  in  the  tapwater  pressure  to  the  auxil- 
iary condenser  which  caused  pressure  changes  of  several 
millimeters     Df      mercury      in   the      s/stem.  To      avoid      this 

problem,  the  flow  through  the  auxiliary  condenser  had  to  be 
continuously  monitored,  unless  data  were  being  taken  la+e  at 
night  when  there  was  no  demand  on  the  laboratory  building's 
water  supply.  The  test  tube  could  be  easily  monitored 
through  the  viewport  for  confirmation  of  filmwise  condi- 
tions. If  there  was  any  sizeable  change  to  dropwise 
condensation  on  the  tube,  the  data  set  was  discarded  and  the 
procedure   was   repeated. 
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III.    FILMS ISE    CONDENSATION 

A-       THE    DROPWISE    PROBLEM 

It  was  essential  during  the  course  of  this  thesis  to 
collect  data  under  filmwise  conditions.  Numerous  problems 
were  encountered  by  Graber  (Ref.  7]  in  avoiding  the  tran- 
sition zo  dropwise  condensation  during  operation,  and  his 
proposed  solution  was  a  vacuum-tight  test  apparatus.  Since 
the  tube  surface  would  wet  completely  after  installation, 
contaminants  leaking  into  the  system  were  possibly  adhering 
to   the   surface    and    promoting   the    dropwise   condensation. 

However,  even  after  obtaining  a  vacuum-tight  apparatus, 
the  author  was  still  unable  to  maintain  good  filmwise 
condensation  for  more  than  two  hours  on  a  smooth  tube,  while 
this  was  enough  time  to  collect  a  complete  set  of  data  for 
this  tube,  filmwise  condensation  lasted  seventes.i  minutes  at 
the  most  on  any  of  the  finned  tubes,  the  average  time  being 
less  than  ten  minutes.  This  was  predictable  -  the  corners  of 
the  fin/surface  interfaces  provided  a  better  trap  for 
contaminants     and   were   harder   to   cl2an   -    but    unacceptable. 

The  use  of  the  steam-cleaning  method  described  in 
Appendix  C  would  thoroughly  clean  the  tube  so  that  complete 
filmwise  condensation  was  re-established,  but  dropwise 
condensation  would    again  become   prevalent   within    minutes. 

B.       SOLUTION 

Having  eliminated  the  possibilities  of  iistalling  a 
dirty  tube  or  contamination  due  to  leakage,  the  only  reason 
for  the  dropwise  problem  had  to  be  coming  from  oitgassing  of 
the  nylon  holders  for  the  test  tube.  The  outgassing  rate  for 
nylon   was   found    to    be  almost   two    oriers   of   magnitude   greater 
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than  the  rate  for  Teflon  [Ref.  9],  Teflon  beia  g  the  same 
order  as  stainless  steel.  By  punping  down  the  apparatus 
while  it  was  hot,  nylon  molecules  were  being  outgassed  into 
the  test  section  and  were  immediately  being  deposited  onto 
the  surface  of  the  cooler  test  tubs,  resulting  in  inadver- 
tant  "sput-ering"  of  the  tube  with  a  nylon  coating  and 
subsequent  dropwise  condensation.  To  eliminate  the  nylon 
interface  with  the  interior  of  ths  tast  section,  special 
stainless  stael  caps  were  manufactured  to  slip  over  the 
nylon   holders. 

Teflon  bushings  were  fitted  within  the  caps  to  insulate 
them  from  the  the  tube.  Teflon  had  both  a  low  thermal 
conductivity  and  a  relatively  low  Dutgassing  rate.  Figure 
3.1    shows   a   detailed   sketch    of    the   installation. 

This  conf igurat icn  appeared  to  solve  the  problem  of 
dropwise  condensation.  Although  an  actual  endurance  test  was 
not  conducted,  the  system  was  in  operation  intermittantly 
for  over  fifteen  hours  with  the  smooth  tube  and  over  four 
hours  with  each  finned  tube  with  no  breakup  of  the  filmwise 
condensation. 

Once  installed,  this  arrangemant  also  eliminated  the 
need  for  the  tube-cleaning  procedure  recommended  by  Graber 
[Ref.  7],  who  felt  that  a  strong  claaning  solution  of  sodium 
hydroxide  and  ethancl  was  needai  to  decontaminate  the 
surface.  Only  a  warm  solution  of  Sparkleen  was  used 
throughout    the    data-ccllection    staga    of   this   thesis. 
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C.       VAPOB    7EL0CITI    LIMITATIONS 

To  check  "che  outside  haat- transf ar  coefficia.it  data  with 
ths  Nusselt  prediction,  and  also  to  obtain  an  accurate 
inside  coefficient  with  the  modifiad  Wilson  Plot,  vapor 
velocities      approaching   zero      were   preferred.  Due    to      the 

design  of  the  system,  however,  the  pressurs  drop  to  the 
auxiliary  condenser  required  vapor  flow  past  the  test  tube. 
This  being  the  case,  attempts  wers  made  to  rainimizs  this 
flow  velocity  by  cutting  down  the  power  to  -hs  boiler  and 
throttling  back  on  the  cooling  water  supply  to  ths 
condenser.  As  vapor  velocity  was  decreased,  however,  drop- 
wise  condensation  again  took  place  on  the  test  tube.  Under 
atmospheric  pressure,  this  occurrence  was  at  vapor  veloci- 
ties of  about  0.5  m/s  (1.6  ft/sec),  and  under  vacuum 
operations  it  occurred  at  about  0.9  m/s  (3.0  ft/sec). 
Apparently  thsre  still  existed  a  sizeable  rate  of  outgassing 
within  the  test  section,  those  gases  collecting  about  the 
test  tube  and  interfering  with  the  filmwise  condensation 
process.  This  suspicion  was  confirmed  when  an  increase  in 
the  vapor  velocity  eliminated  the  formation  of  drops  on  the 
tube . 
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IV.    DATA    MMCriON 

The  data  collected  and  stored  on  the  computer  disks 
reduced  using  the  programs  WILSON,  SIEDER,  and  DRP.  The 
programs  were  amenable  to  changes,  which  allowed  the  author 
to  analyze  and  compare  results  while  varying  parameters 
within  the  programs.  Stepwise  reduction  procedures  for  each 
program  are  listed  below  and  the  program  listings  are  found 
in   Appendix   E. 

A.  FBOGBAM    SIEDEB 

1.  Compute    the   average    cooling    water   temperature. 

2.  Compute    the    average    wall   tenperature. 

3.  Compute    -che    cooling    water   valocity. 

4.  Compute   the    mass-flow    rate   of    the   cooling    water. 

5.  Compute    the    heat   transferrai    to    the   cooling    water. 

6.  Compute   the   average    inside   wall    temperature. 

7.  Compute   the    log-mean-temperature    difference. 

8.  Compute    the    inside    heat- transfer    coefficient. 

9.  Compute    the    Nusselt    number. 

10.  Compute    the    Sieder-Tate   parameter. 

11.  Compute    the    inside    coefficient. 

B.  PBOGSAH    WILSON 

1.  Assume    a   value  for    the   Sieder-Tate   coefficient. 

2.  Compute  the  Reynolds  and  Prandtl  numbers  for  flow 
through    the   tube. 

3.  Compute  the  log-mean-temperature  difference,  heat 
flux,      and   overall    heat-transfer    coefficient    for    the    tube. 

4.  Assume   an  outer  tube   surfacs    temperature. 
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5.  Compute      the   outside      condaasing      coefficient      using 
properties    evaluated   at   the    film   tsiperature. 

6.  Compute      the   outer      surface   temperature      and    iterate 
steps   5  and  6    if   not   within    1%. 

7.  Assume        a      viscosity-correction        factor      for        the 
Sieder-Tate   equation    cf    1.0. 

8.  Compute    the    inside    heat-transfer    coefficient- 

9.  Compute   the    inner   surface   tsmperature. 

10.  Compute    the    viscosity      correction   factor    and   iterate 
steps   8  through    10    if  not   within    ^%. 

11.  Compute        the  Sieder-Tate        coefficient    and      iterate 
steps    2  through    11    if  not   within   0.5%. 

C.  PROGRAH    DEP 

1.  Compute    the    average    ceding    wazer    temperature. 

2-  Compute   the    cooling    water   vslocity. 

3.  Compute    the    mass-flow   rate   of    the   cooling    water. 

U.  Compute    the    heat   transferred    to    the   water. 

5.  Compute    the    log-mea n-tempsrat ure    difference. 

6.  Compute    the    overall    heat-"!:ransf er    coefficient. 

7.  Compute    the    wall  resistancs    of    the   tube. 

8.  Compute   the    Reynolds   number   of   the   cooling   water. 

9.  Compute   the    inside    heat-transfer   coefficient. 

10.  Compute    the    condensing   heat-transfer    coefficient. 
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^«     RESULTS    AND    DISCUSSION 

Numerous  data  runs  wera  made  using  ths  procedure 
described  in  Section  II. G.  Time  constraints,  however, 
limited  the  number  cf  repeat  runs  than  could  be  made  for 
this  thesis.  Primary  concern  focused  on  establishing  a  reli- 
able, repeatable  Sieder-Tate  coefficient.  Daca  were  taken 
for  all  of  the  -cubes  at  both  atmospheric  and  vacuum  (88mni 
Hg,  1.7  psia)  conditions.  Complete  filmwise  condensation  was 
maintained  for  all  data  runs,  and  the  mass  concentration  of 
noncondensable  gases  was  held  between  0%  and  -1%  during  all 
testing.  The  negative  value  was  indica-cive  of  slight  super- 
hear    in    rhe   system   or  an   inaccurate   manometer   reading. 

A.       INSIDE    HEAT    TRANSFER    COEFFICIENT 

"^  •      InsT-rumented    Tube   R  esults 

Figure  5.1  shews  the  variation  of  the  Nusselt  number 
as  a  function  of  the  Sieder-Tate  parameter  for  the  instru- 
mented -cube  run  at  atmospheric  pressure  (run  SDA7)  .  This 
method  Yi®l<i®d  ^  Sieder-Tate  coefficient  of  0.035  on  two 
separate  data  runs  (SDA7  and  SDA8)  .  The  same  method  under 
vacuum  conditions  yielded  a  coefficient  of  0.037,  which  was 
also  repeated  (runs  SDA5  and  SDA6) .  The  temperature  distri- 
bution around  the  tube  wall  was  symmetrical  about  the 
vertical  plane  passing  downward  through  the  centerline  of 
the  tube,  and  showed  as  much  as  a  15  C  temperature  drop 
from    the   top  of    the    tube   to    the   bottom. 

2-      Smooth   Tube    Results 

Figure  5.2  shows  the  modified  Wilson  Plot  for 
smooth-tube      data     collected      at    atmospheric      pressure      (run 
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0D12).  This  method  produced  a  Sisdar-Tate  corfficien*.  of 
0.033  with  an  intercept  of  1.16.  A  similar  plot  for  vacuum 
conditions  (run  OD10)  produced  a  coefficient  of  0.036  and  an 
intercept  of  1.27.  The  Nusselt  theory  predicts  an  intercept 
of  1.53  but,  as  was  noted  in  Section  III.C,  a  finite  vapor 
velocity   past    the  tube   was    unavoidable   in   this  research. 

Figure  5.3  shows  a  plot  of  rhe  Fujii  [ Hef .  10] 
correlation  for  the  smooth-tube  rans.  This  experimental 
correlation  is  to  correct  the  Nusselt  number  for  the  effects 
of  vapor  shear  on  the  test  tube  and  is  a  plot  of  the  equa- 
tion : 

NuRe-i/2    =   0.96FI/5 

The  results  obtained  were  slightly  higher  than  those 
predicted  using    this    correlation. 

3  •      Discrepa  ncie  s 

The  original  Sisder-Tate  aquation  [Raf.  11]  for 
fully  developed  turbulent  flow  in  a  rube  with  an  L/D  ratio 
greater  than  60  has  a  leading  cosfficient  of  0.027,  so  a 
higher  value  for  the  shorter  tube  (L/D  =  18)  used  for  this 
research  stands  tD  reason.  The  data  for  the  tjbes  appears 
at  first  inconsistent,  with  results  of  0.033,  0.035,  0.036, 
and  0.037  -  showing  up  to  a  6%  scatter  from  the  aean.  Both 
tubes,  however,  showed  a  largar  coefficient  when  the  test 
section  pressure  was  reduced,  which  reduced  the  heat  flux  as 
well.  Dropping  the  pressure  from  atmospheric  to  a  high 
vacuum  (88  ram  Hg,  1.7  psia)  reduced  the  heat  flux  by  a 
factor  of  almost  three  [Fig.  5.4].  Stated  another  way,  a 
low  heat  flux  produced  coefficients  of  0.036  and  0.037, 
while   a    higher    heat    flux   produced   0.033    and   0.035. 

The  author  doubts  the  reliability  of  the  inside  co- 
efficient obtained      under  vacuum   conditions    for   two    reasons: 
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First,  reducing  the  flow  rate  through  the  tube 
lowers  the  heat  flux  and  increases  the  inside  ccnvective 
thermal    resistance.  Since   the      condensing    coefficient      is 

inversely  proportional  to  the  heat  flux  raised  to  the  one- 
third  power,  the  outside  condensing  resistance  decreases. 
Because  of  this,  the  inside  resistance  becomes  the  dominant 
factor  in  the  overall  thermal  resistance  of  the  tube.  Any 
error  m  measuring  the  inside  hea-c  transfer  coefficient  will 
greatly  amplify  the  computed  error  in  the  condensing  resis- 
tance. Thus,  the  condensing  coefficient,  is  highly  sensitive 
to  the  accuracy  of  the  inside  coefficient,  particularly  at 
low    water   velocities    (and  corresponding    heat    fluxes)  . 

Second,  the  test  apparatus  was  subject  to  radio- 
frequency  (RF)  interference  from  the  power  supply.  As  the 
power  was  cut  back  to  lower  the  heat  flux,  the  silicon 
controlled  rectifier  chopped  the  input  signal  proportionally 
and      emitted  the      chopped   portion      as      RF    enargy.  Despite 

efforts  to  shield  the  cabling  to  the  data- acquisition 
system,  it  still  captured  the  HF  signal  with  unacceptable 
results,  and  the  thermopile  was  experiencing  instantaneous 
fluctuations  of  several  hundred  microvolts.  Increasing  the 
power  supply,  however,  decreased  ths  errant  signal  until,  at 
the  maximum  power  input  used  for  th=  high  heat  flux  measure- 
ments, the  thermopile  showed  a  near-steady  temperature 
difference   that    was    within    0.1    °C   of    the   quartz   thermometer. 

For  these  reasons,  the  Sieder-Tate  coefficients 
obtained  from  the  low  heat-flux  runs  were  discarded. 
Averaging  the  coefficients  found  from  both  measurements  for 
the  high  heat  flux  case  yielded  a  coefficient  of  0.034  ± 
0.001. 
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B.        FINNED    TUBES 

Figure  5.5  shows  a  relativa  plot  of  the  outsida 
condensing  coefficient  for  all  six  finned  tubes  compared 
with  smooth  tube  results  and  the  Nussalt  line.  Figure  5.6 
shows  the  same  results  corrected  for  area  increases  due  to 
finning.  These  data  were  computed  using  the  leading  co<=iffi- 
cient  on  the  inside  of  0.034  obtained  in  the  previous 
sec- ion. 

All  tubes  tested  show  good  enhancement  over  the  sniooth 
tube.  Notice  the  increase  of  the  smooth  tube  over  the 
Nusselt  line  -  a  resulx.  of  the  inherent  vapor  shear  in  the 
test  section.  The  increased  scat-er  in  the  lower  heat-flux 
range  is  agaii  a  function  of  the  increased  significance  of 
errors  in  measuring  the  inside  thermal  resistances.  The 
plotting  program  used  neglected  ordinate  values  below  zero 
or  greater  than  10^,  so  the  scattering  in  this  portion  of 
the  plot  is  actually  worse  than  it  appears.  Figure  5.7  is  a 
plot  of  the  values  obtained  for  the  outside  heat  transfer 
coefficients  of  the  six  finned  tubes  for  a  constant  heat 
flux  of  250,000  W/m2  (79,000  Btu/hr-ft2)  and  a  pressure  of 
88  mm  Hg  (1.7  psia) .  This  plot  more  clearly  shows  the 
optimum    pitch   of    2.5    mm. 

The  heat-transfer  characteristics  for  any  tube  will  be 
enhanced   by     the  addition  of   fins.  In    the    case      of   purely 

convective  heat  transfer,  the  enhancement  is  a  function  of 
the    increased  surface  area    exposed   to    the    fluid   medium. 

Filmwise  condensation,  however,  is  dissimilar  in  that 
the  build-up  of  a  condensate  fila  acts  as  an  additional 
thermal  resistance  for  the  heat-transfer  process.  The  objec- 
tive of  tube  enhancement,  therefore,  is  to  decrease  this 
film  thickness.  The  surface  tension  forces  in  the  conden- 
sate tend  to  draw  the  liquid  toward  the  fins,  leaving  the 
tube  surface  with  a  thinner  film.  The  thinner  film  results 
in   a    higher   haat-transfer   coefficient.      [Bef.    12]. 
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The  optimum  pitch  found  in  this  research  increases  the 
surface  of  the  condenser  tube  88'^,  but  condensing  heat- 
transfer  coefficients  for  that  pitch  were  enhanced  by  as 
much    as    330%. 

The  optimum  pitch  serves  as  the  right  trade-off  between 
the  attractive  forces  of  the  fins  on  the  condensate  and  the 
channel  area  between  the  fins  to  drain  the  condensate  from 
the  tube.  A  pitch  smaller  than  th=  optimum  has  too  narrow  a 
gap  between  fins  to  efficiently  allow  condensate  run-off. 
This  is  because  the  condensate  drawn  to  adjacent  fins  only 
combines   to   create    a   thicker    film   as    shown    in   Figure    5.8. 
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Figure  5.8       Sketch   of  the    effect  of    finning   on   condensation, 
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VI.     CONCLUSIONS    AND   SECOMMENPATIONS 

A.  CONCIOSIONS 

1.  The  test  apparatus  successfully  operates  under 
vacuum  conditions  with  no  degrading  effects  caused  by  the 
presence   of  nDncondensable    gases. 

2.  The  Sieder-Tate  coefficient  was  deremined  to  be 
0.034  ±  0.001  using  both  an  instrumented  smooth  tube  and  the 
modified   Wilson    plot    technique. 

3.  An  op-cimum  fin  pitch  of  2.5  mm  was  found  for  the 
finned   tubes  tested. 

U.  Heat-transfer  enhancements  up  to  3007,  were  raalized 
while  the  increase  in  the  tube  outer  surface  area  (due  to 
finning)    was  only  88%  over    the    smooth   tube. 

B.  HECOMHENDATIONS 

1.  Install  a  second,  more  reliable  pressure  transducer 
to   replace   the    manometer. 

2.  Install  a  larger  cooling-wat sr  sump  to  eliminate  use 
of  the   once-through    system    from   the    frash-water   tap. 

3.  Collect  data  using  an  insert  for  the  test  tube  to 
decrease  the  internal  convective  resistance  and,  thereby, 
decrease  the  uncertainty  m  the  calculated  condensing  coef- 
ficient. 

U.  Take  data  for  the  finned  tubes  while  varying  fin 
height   and   thickness    (as   well   as   pitch). 
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APPENDIX    A 
THERMOCOOPLE    CALIBRATION 


1.  EQUIPMENT    USED 

a.  Thermocouple    Wire 

Copper-constantan,  0.2U5-ra!n  (0.01-in)  Teflon- coated  wire 
was    used   for  all  theriDOcouples. 

b.  Calibration    bath 

A  Rosemont  Engineering  Model  91 3A  calibration  bath  was 
used.  A  schematic  representation  of  the  bath  is  shown  in 
Figure   A.I. 

1)  Heating:      Elecrrical 

2)  Cooling:      Liquid   Nitrogen 

Note:  Once  a  desired  temperature  is  reached,  the  temp- 
erature is  held  constant  by  rapid  cycling  between  heating 
and  cooling.  The  bath  is  rated  for  temperature  fluctuations 
of   less   than  O.OI'^C. 

c.  Thermocouple    readout 

An  HP  3054A  Data  Acquisition/Control  System  was  used  to 
obtain   data.      Resolution  of    the   acquisition   system    was    I^V. 

d.  Eath  temperature   measurement 

A  platinum  resistance  thermometer  with  an  accuracy  of 
0.01      C    was   used. 

2.  PREPARATION 

a.  Procedure 

The  instrumented  tube  (with  the  wall  thermocouples 
installed)  and  the  steam  thermocouples  were  immersed  in  the 
bath    as    well  as    the    probes    for   the   quartz   thermometer. 

b.  Analysis 

The  computer  program  TCAL  was  used  to  monitor  and  store 
all  thermocouple  readings  on  a  disk.  A  listing  of  the 
program   is   located    in  Appendix   E. 
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3.       CALIBRATION    PROCEDURE 

a.  The   bath   temperature   was   set    at    about    10 ''C. 

b.  When  the  bath  temperature  reached  steady  s-^-ate,  its 
value    was   entered  into  the    computer. 

c.  The  computer  automatically  recorded  and  printed  all 
thermocouple  readings. 

d.  The  bath  temperature  was  raised  in  incrsments  of  10 
C  to  90  °C  and  steps  b  and  c  were  repeated  for  each  incre- 
ment . 

a.       CALIBRATION    CURVES 

a.  A  least-squares  method  was  used  to  generate  a  poly- 
nomial  of   the   form: 

Dr   =    a^    *   a,.T      +    a^TS 

where:  D^  is  the  difference  betwsen  the  bath  and  thermo- 
couple  temperatures,    and 

T  is  the  value  of  the  thermocouples  obtained  using  the 
seventh-order  polynomial  fitted  for  the  Type  T  thermocouple 
wire  used.  This  pclynomial  is  listed  in  the  beginning  of 
the    program  TCAL. 

b.  Coefficient    values    obtained    were: 

a^    =    U.7338  x    10-3 
a,    =    7.6928   x    10-3 

a^  =    -8.0779  X    10-5 

c.  A  plot  of  the  curves  is  found  in  Figure  A.I.  The 
curve  for  the  wall  thermocouples  gives  a  different  reading 
due  to  the  thermal  conduction  of  ths  air  temperature  through 
the  tube.  The  data  used  for  calibration  did  not  include 
this  average  value.  Only  the  value  obtained  for  the  two 
steam   thermocouples    was   used   for   all    of    the   thermocouples. 
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APPENDIX    B 
SAMPLE    CALCULATIONS 


1.        HEAT    LOSS 

A  set  of  calculations  was  parfornisd  to  estimate  the  heat 
loss  of  the  apparatus  due  to  natural  convection.  These 
calculations  were  performed  for  two  main  sections:  the 
boiler  (with  and  without  insulationi  and  the  piping  to  the 
test  section.  The  minimum  power  input  used  during  the  thesis 
was    17    ;^W. 

a.      Boiler    (without   insulation) 

1)  Given      the    following      dimensions    and      properties   of      the 
toiler  section: 

r,',5   =    0.152U   m 

^c*  3    =    0. 1 56 1    m 

L  J       =    0.406a    m 

k  (^       =1.4    W/m-K 

2)  Given  the   following   temperatures: 
T„B   =    lOO^C 

Tc..s=  95°C 
T.^  =  25*C 
T^      =    (95+25)/2    =   60^C 

3)  Given   the   following      properties   of   air   at   the    film   temp- 
erature: 


k^     =   28.74x10-3    W/m-K 
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Pr      =    0.702 


-t^      =    19.21x10-*      mz/s 


0       =    0.003    K-i 

4)  Compute     ths   Srashof        Number   for        the   boilsr      assuming 
laminar  conditions    (Ra  <   10^): 

Grg  =   g(3(T,  g  -T^)L3/^2 
Grg    =    3.7x108 

5)  Compute   the    Rayleigh  Number   for   rha    boiler: 

Rag   =  Gr  Pr 

Ra^j  =   2.6x108  (<    109) 

6)  Determine        the    applicability      of   a      flat-plate    analysis 
according  to  the   method   of    Sparrow   and   Gregg   (Ref-    13] 

L  /r„,e  <      (0.15/>r5)Gr»/* 

2.60        <      7.36 

8)  Compute  the    average    Nusselt   Numbar    (via   flat-plate   anal- 
ysis) : 

Nu      =    EO.825    ♦    0.387Rag  i/*/[  1  +  (a.492/?r) 9/1* ]8/27]2 
Nu      =   81.3 

9)  Compute  the    external  heat-transfer   coefficient: 
ho    =    Nu   k   /L 

he,    =    5.75    W/m2-K 

10)  Compute  the    external   convective   thermal  resistance: 
R      =    1/(h^2Trro,.3L^) 

R      =0.436    K/W 

11)  Calculate  the  wall  resistance: 
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Rc«  =   In  (r^^B/r,,g) /2-7tJc^Lb 

R^v  =   0.007    K/W 

12)     Neglecting    the    wall   resistance, 
from    the   boiler: 

Qg  =    (T.>-T^)/(R    ♦Ew) 
Q  0  =    169.3    W 

b.  Boiler    (with   insulation) 
Qb  =    29.2    W 

c.  Piping 

Q  p  =  65.0   W 
a.      Total    (worst   case) 


calculate  -^he    heat   less 


Q      =   Q  3  >   Q  P 
Q      =    234.3    W 


(<<   17000   W    of    input   power) 
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APPENDIX    C 
STEAM-CLEANING    PR3CED0HE 

Note:  The  procedures  listed  hare  provide  an  excellent 
method  of  cleaning  the  apparatus,  but  will  cause  the  walls 
of  the  apparatus  to  heat  to  temperatures  over  100  C  (212  F)  . 
Once  this  happens,  several  hours  are  required  before  th s 
walls  can  cool  down  such  that  they  won't  superheat  the  staam 
generated  during  an  experimental  run.  Therefore,  use  this 
method  only  prior  to  operating  at  atmospheric  pressure  or 
when    the  system    is    badly  contaminated. 

1.  Ensure  the  water  level  in  the  boiler  is  at  leasT  six 
inches  above  the  upper  ends   of   the   heaters. 

2.  Remove  the  thermopile  from  the  inlet  side  of  the 
test    section, 

3.  Energize  the  power  supply  and  adjust  the  rheostat 
until   the   voltmeter   reads  about   100    V. 

U.      Close  the  recirculating  valve. 
5.      Open   the   fill/drain    valve. 

6-  Once  vapor  begins  exiting  the  fill/drain  line,  allow 
the    sys-em   to   steam    for   several   minutes. 

7.  To  maximize  the  steam  flow  through  the  test  section, 
fully    open    the    throttle   to    the   auxiliary   condenser. 

8.  When  done,  simultaneously  close  the  drain/fill  valve 
and    open  the  recirculating    valve. 

9.  Circulate  cooling  water  through  the  test  tube  and 
check   for  the  presence  of  any    dropwise   condensation. 

10.      Reinstall  the  thermopile. 
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APPENDIX    D 
SYSTEM    STAET-OP    AND    SHOT-DOWH    PROCEDURES 

To   start   the  system: 

1.  Fill  the  boiler  to  at  least  six  inches  above  the 
upper    level  of    the    heaters, 

2.  Energize  the    air   ejector   for    ten   minutes. 

3.  Circulate  water  through  tha  cooling  water  sump  by 
opening   -he   inlet   valve   and    activating   the    siphon. 

1.  Energize  the  power  to  the  boiler  and  adjust  the 
rheostat  until  the  voltmeter  reads  90  V  for  vacuum  opera- 
tion,  or    170  V    for    atmospheric    pressures. 

5.  Energize  the  circulating  pumps  and  adjust  the 
throttle   for  the   desired   flow   through   the   tube. 

6.  Open  the  valve  to  the  auxiliary  condenser  and  adjust 
the  flow:  about  15%  for  vacuum  operation  or  30%  for  atmos- 
pheric runs. 

7.  Observe  the  steam  temperature  indicated  on  th  efront 
of  the  data  acquisition  unit  until  the  voltage  is  reached 
corresponding  to  the  saturation  temperature  Df  the  the 
desired  operating  pressure. 

8.  Energize  the   air   ejector   for    two    minutes. 

9.  Adjust  the  the  flow  rate  through  the  auxiliary 
condenser  to  bring  the  apparatus  to  the  desired  operating 
pressure. 

To   secure  the  system  : 

1.  Secure   cooling  water  to   the   test   tube. 

2.  Secure    power    to   the    boiler. 

3.  Open  the   vent  valve. 

U.      Secure   the    water   supply   to  the    auxiliary    condenser. 
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APPENDIX    E 
COMPUTER    PBOGBAH    LISTINGS 

The   following      pages   contain   the    computer      program    list- 
ings   used   for   this    thesis: 

SIEDER  (page   61) 

WILSCN  (page  65) 

DRP  (page  7  1) 

TCAL  (page   80) 
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PROGRAM    SIEDER 


1000!  FILE  MAME:    SIEDER 

1010?  DISK  NUMBER:  12 

1020?  REVISED:     November  29.  '983 

1030? 

1040  COM  /Cc/  C(7) 

1050  DIM  EmfdO)  .Tuj(5) 

i060  DATA  0. 10086091  .25727.94369.-757345. 8295. 78025595. 81 

1070  DATA  -9247486589. 6. 97688En .-2. 56192E13. 3. 94078E14 

T080  READ  C(*) 

1090  Kcti=385 

1100  Di=.0127 

1110  Do=. 01905 

1120  Dr=. 015875 

1130  L=. 13335 

1140  L1=. 060325 

1150  L2=. 034925 

1160  PRINTER  IS  701 

1170  BEEP 

1180  CLEAR  709 

1190  INPUT  "ENTER  MONTH,  DATE  AND  TIME  (MM: DD:HH:MM :SS)" .BS 

1200  OUTPUT  709:"TD":BS 

1210  Series:? 

1220  OUTPUT  709:"TD" 

1230  ENTER  709:AS 

1240  PRINT  USING  "1  OX .""Month,  date  and  time:  "".14A":AS 

1250  BEEP 

1260  INPUT  "ENTER  DISK  NUMBER", Dn 

1270  PRINT 

1280  PRINT  USING  "10X,""NOTE:  Program  name:  SIEDER 

1290  PRINT  USING  "16X.""Disk  number   =  "".DD";Dn 

1300  BEEP 

1310  INPUT  "ENTER  INPUT  MODE  (  1  =305^A  .2=FILE)" , Im 

1320  IF  Im=1  THEN 

1330  SEEP 

1340  INPUT  "GIVE  A  NAME  FOR  THE  DATA' FILE" .D_fi leS 

1350  CREATE  BDAT  D_fileS.10 

1360  ELSE 

1370  BEEP 

1380  INPUT  "GIVE  THE  NAME  OF  THE  DATA  FILE" .D_f i leS 

1390  BEE? 

1400  INPUT  "ENTER  THE  NUMBER  OF  RUNS  STORED". Nrun 

1410  PRINT  USING  "16X,""This  analysis  uas  performed  for  fxle  "" .  1  OA" ; D_f i Ie$ 

1420  END  IF 

1430  BEEP 

1440  INPUT  "GIVE  A  NAME  FOR  PLOT  DATA  FILE". Plots 

1450  BEEP 

1460  INPUT  "ENTER  OPTION  FOR  END-FIN  EFFECT  (1 =Y , 0=N)" . Ife 

1470  IF  lfe=0  THEN  PRINT  USING  "16X.""This  analysis  neglects  end-fin  effect"' 

1480  IF  Ife=l  THEN  PRINT  USING  "16X.""This  analysis  includes  end-fm  effect"' 

1490  CREATE  BDAT  Plots. 5 

1500  ASSIGN  »File  TO  D_fileS 

1510  ASSIGN  §Filep  TO  Plots 

1520  J=0 

1530  Sx=0 

1540  Sy=0 

1550  Sxs=0 

1560  Sxy  =  0 

1570  IF  Im=1  THEM 

1580?  READ  DATA  THROUGH  THE  DATA  ACQUISITION  SYSTEM 
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1590!  IF  THE  INPUT  MODE  (Im)  =  1 

1600  BEEP 

1610  INPUT  "ENTER  FLOWMETER  READING", Fm 

1620  OUTPUT  709;"ftR  AF20  AL30  \/R1" 

1630  FOR  1=0  TO  10 

1640  OUTPUT  709:"fiS  SA" 

1650  IF  I>4  THEN 

1660  Se=0 

1670  FOR  K=0  TO  19 

1680  ENTER  709:E 

1690  Se=Se+E 

1700  NEXT  K 

1710  Emf<I)=ABS(Se/20) 

1720  ELSE 

1730  ENTER  709 ;E 

1740  Emf<I)=ABS(E) 

1750  END  IF 

1760  NEXT  I 

1770  OUTPUT  713:"T1R2E" 

,1780  WAIT  2 

1790  ENTER  7^3^,J]] 

1800  OUTPUT  713;"T2R2E" 

1810  WAIT  2 

1820  ENTER  713:T2 

1830  OUTPUT  713;"T1R2E" 

1840  WAIT  2 

1850  ENTER  713:T12 

1860  T1=(T11+T12)».5 

1870  ELSE 

1880?  READ  DATA  FROM  A  USER-SPECIFIED  FILE  IF 

1890?  INPUT  MODE  (Im)  =  2 

1900  ENTER  «File;E!nf(»).T1  .T2.Fm 

1910  END  IF 

1920  Tavq=<T1+T2)-.5 

1930  Tuiail=0 

1940  FOR  1=5  TO  10  ^ 

1950  Tu<I-5)=FNTvsv(Emf(I)) 

1970  T(jail  =  Tuiall  +  Tuj<I-5) 

1980  NEXT  I 

1990  Tuiall  =  Tu)ail/6 

2000  Cpu=FNCpu(Tavg> 

2010  Rhou=FNRhoij(Tavg) 

2020  Md=5.00049E-3+6.9861937E-3»Fin 

2030  Md  =  Md»(1  .0365-1  .96644E-3-T1 +5  .252E-6*Tr2)/ .995434 

2040  Mf=Md/Rhow 

2050  Vu=Mf/(PI»Di^2/4) 

2051  T2c  =  T2-(.0138+.001»Vu)~2) 
2060!  T2c=T2-.004»yu'2 

2070  Q=Md»Cpu*<T2c-T1) 

2080  Dtuj=Q»L0G<Do/Di)/(2*PI»Kcii»L)*.5 

2090  Tuall=Tuall-Dtu 

2100  Lmtd=(T2c-T1)/LQG(<Twail-T1)/(TuaiI-T2c)) 

2110  Ku=FNK(j(Tavg) 

2120  P1=PI»(Di+Do) 

2130  P2=PI»<Di+Dr> 

2140  A1=<Do-Di)»PI«Do 

2150  A2=(Dr-Di)»PI»Dr 

2160  Hi=Q/(PI»Di'L»Lmtd) 

2170  IF  lfe=0  THEN 

2180  Hic=Hi 

2190  GOTO  2300 
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2200 
22!0 
2220 
2230 
22^0 
2250 
22S0 
2270 
2280 
2290 
2300 
2310 


END  IF 

M1=(Hi»P1/<Kca*AI ))  '.5 

M2=(Hi-P2/(Kca*A2))".5 

Fel =FMTanh(M1»LI )/<M1»Ll ) 

Fe2=FNTanh<M2»L2)/(M2*L2) 

Hic=Q/<PI»Di*(L+LI*Fel+L2»Fe2>»Lmtd) 

IF  ABS(<Hi-Hic>/Hic)>.01  THEN 

Hi=(Hic+Hi)».5 

GOTO  2210 

END  IF 

PRINT 

PRINT  USING  "1  OX, ""Position  number 
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2320  PRINT  USING  "10X.""Vlail  tenperatare  (Dea  C)  :  ""  .6(DD  .  DO  ,  1 X)"  :  Tu»(*) 

2330?  CALCULATE  THE  NUSSELT  NUMBER 

2340  Nu=Hic»Di/Ku 

2350  Mau=FNMiiu<Tavg) 

2360  Re=Rhou»^;u)*Di/f1atj 

2370  Cf  =  <Muw/FNf1uu(Tuall>)' .14 

"2380  Pru)=FNPru»<Tavg) 

2390  X=Re^.8»Pru-.3333»Cf 

2400?  COMPUTE  COEFFICIENTS  FOR  THE  LEAST-SOUARES-FIT 

2410?  STRAIGHT  LINE 

2420  OUTPUT  SFilepiX.Na 

2430  PRINT  USING  "10X.""Tuall    Tin    Tout   Lmtd   (J<j  X        Nu 

2440  PRINT  USING  "1  OX .4(20 .2D .2X) ,Z.DD,2X .50 .D ,2X .4D .0" : Tuai 1 ,T1 , T2c .Lmtd . Vw.X . 

Nu 

2450  Sx=Sx+X 

2460  Sy=Sy+Nu 

2470  Sxs=Sxs+X»X 

2480  Sxy=Sxy+X»Nu 

2490?  STORE  RAW  DATA  IN  A  USER-SPECIFIED  FILE  IF 

2500?  INPUT  MODE  dm)  =  1 

2510  IF  Ini=1  THEN  OUTPUT  ^Fi  le:Emf  <•)  ,  T1  .  T2.Fm 

2520  BEEP 

2530  J=J+1 

2540  IF  Im=1  THEN 

2550  INPUT  "ARE  YOU  TAKING  MORE  DATA  ( 1 =YES.0=N0)7".Go_on 

2560  Nrun=J 

2570  IF  Go_on=;  THEN  1570 

2580  ELSE 

2590  IF  J<Nrun  THEN  1570 

2600  END  IF 

2610?  Ci=Sxy/Sxs 

2620  Ci=<Nrun»Sxy-Sy*Sx)/<Nrun*Sxs-Sx  2) 

2630  Ac=(Sy-Ci»Sx)/Nrun 

2640  PRINT 

2650  PRINT  USING  "10X. ""Sieder-Tate  Coefficient   =  "",D.4D";Ci 

2650  PRINT 

2570  PRINT  USING  "1  OX .""Least-Squares  Line: 

2680  PRINT  USING  "12X. ""Slope        =  "" ,MD.5DE ." :Ci 

2690  PRINT  USING  "12X , ""Intercept    =  ""  ,MD.5DE, " :Ac 

2691  PRINT 

2700  IF  Im=1  THEN 

2720  BEEP 

2730  PRINT  USING  "10X,""NOTE:  "",11.""   data  runs  were  stored  in  file  "",8A";Nru 
n.O_f lies 

2731  ELSE 

2732  PRINT  USING  "10X.""NOTE:  The  above  analysis  was  performed  for  file  "".t4A" 
;D_f ileS 
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2740  END  IF 

2750  PRINT  USING  "MSX .  •'"Plot  data  are  stored  in  file  "•M^A" :  Plots 

2760  ASSIGN  ^File  TO  * 

2770  ASSIGN  aPilep  TQ  * 

2780  END 

27S0  DEF  FNRhou.<T> 

2800  Ro=100G.35724-T*( . 774489-T»(2 .262459E-2-T»3. 03304E-4> ) 

2810  RETURN  Ro 

2820  FNEND 

2830  DEF  FNPru)(T) 

2840  Prtj=FNCpui<T)«FNHtnj(T)/FNK(j<T) 

2850  RETURN  Pru 

28G0  FNEND 

2870  DEF  FNMuu(T> 

2880  A=247.8/(T+I33.15) 

2890  Muu=2.4E-5*10'A 

2900  RETURN  Muu 

2910  FNEND 

2920  DEF  FNKbid) 

2930  Kij=. 572133504477  +  1  .52770121209E-3-T 

2940  RETURN  Ku 

2950  FNEND 

2960  DEF  FNTvsvCEmf) 

2970  COM  /Cc/  C(7) 

2980  S.iin=C(0) 

2990  FOR  1=1  TO  7 

3000  SujTi=SLun+C(I>»Emf' I 

3010  NEXT    I 

3020  RETURN  Stun 

3030  FNEND 

3040  DEF  FNCpu(T) 

3050  Cpu=(4.21 !20858-T»(2.26326E-3-T»<4.42361E-5+T»2.71428E-7)))»1000 

3060  RETURN  Cp.j 

3070  FNEND 

3080  DEF  FNTanh(X) 

3090  P=EXP<X> 

3100  0=EXP<-X) 

3110  Tanh=<P+Q>/(P-Q) 

3120  RETURN  Tanh 

3130  FNEND 
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PROGEAK  WILSON 


1000?  FILE  MAHE:  WILSCM 

1010?  REVISED:    Decemoer  5.  1983 
1020? 

1030  COM  /Cc/  C(7) 

1040  DATA  0. 10086091 .25727. 94369. -7S7345. 3295. 78025595. 81 

1050  DATA  -9247486589. 6. 97688En .-2. 66132E13. 3. 94078E14 

1060  READ  C<») 

1070  DIM  Emf(4) 

1080  L=. 130175 

1090  L1=. 060325 

1100  L2=. 034925 

1110  Do=. 01905 

1120  Di=.0127 

1130  Dr=. 015785 

1140  Kcu-335 

I  150  Rm=Do-LQG(Do/Di)/(2*Kcu) 

1160  PRINTER  IS  701 

1170  BEEP 

•1180  CLEAR  709 

1130  INPUT  "ENTER  MONTH.  DATE,  AND  TIME  (MM:DD:HH:MM:SS",BS 

1200  OUTPUT  709:"TD";BS 

1210  Jp=0 

1220  OUTPUT  709:"TD" 

1230  ENTER  709 -.AS 

1240  PRINT  USING  "lOX. ""Month ,  date  and  time   :  "".14A":AS 

1250  BEEP 

1260  INPUT  "ENTER  DISK  NUMBER", Dn 

1270  PRINT 

1280  PRINT  USING  "10X.""NOTE:  Prograin  name   :  MJrIILSON 

1290  PRINT  USING  "15X,""Disk  number   =  "".DD":Dn 

1300  BEEP 

1310  INPUT  "ENTER  INPUT  MODE  ( 1 =3054A .2=FILE)" . Im 

1320  IF  Im=1  THEN 

1330  BEEP 

1340  INPUT  "GIVE  A  NAME  FOR  THE  DATA  FILE" .D_f i leS 

1350  CREATE  BDAT  D_file$.l0 

1360  ELSE 

1370  BEEP 

1380  INPUT  "GIVE  THE  NAME  OF  THE  DATA  FILE" .D_f i leS 

1390  PRINT  USING  "16X,""Thi3  analysis  is  for  data  in  file  "" , 1 4A" :D_f ileS 

1400  BEEP 

1410  INPUT  "ENTER  THE  NUMBER  OF  RUNS  STORED". Nrun 

1420  END  IF 

1430  BEEP 

1440  INPUT  "GIVE  A  NAME  FOR  PLOT-DATA  FILE". Plots 

1450  BEEP 

1460  INPUT  "ENTER  OPTION  ( 1 =QCT,2=T-PILE.3=AVE>"  ,  Itm 

1470  BEEP 

1480  INPUT  "ENTER  OPTION  FOR  END-FIN  EFFECT  (1 =Y,0=N)". Ife 

1490  IF  Itm=1  THEN  PRINT  USING  "16X,""This  analysis  ases  OCT  readings 

1500  IF  Itm=2  THEN  PRINT  USING  "16X.""This  analysis  oses  T-PILE  readings 

1510  IF  Itm=3  THEN  PRINT  USING  "16X,""This  analysis  ases  average  of  QCT  and  T-P 

ILE  readings' 

1520  IF  Ife=1  THEN  PRINT  USING  "16X,""This  analysis  mclades  end-fin  effect 

1530  IF  lfe=0  THEN  PRINT  USING  "16X,""This  analysis  neglects  end-fin  effect 

1540  CREATE  BDAT  Plots.  10 

1550  ASSIGN  «Filep  TO  Plots 

1560?  Cia=.040 

1570?  Cil=.028 
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1580  Jj=0 

1590  Ci=.03 

1600  J=0 

1610!  Ci  =  <Ci.j.+Cil)».5 

1620  Sx=0 

1530  Sy=0 

1640  Sxs=0 

1550  Sxy=0 

1660  PRINT 

1670  PRINT  USING  "1  OX .""Iteration  number  =  "".DD";Jj  +  l 

1580  IF  Jj=0  OR  Jp=1  THEN 

1690  PRINT 

1700  PRINT  USING   "12X.""T1  T2  Tsat  Lnitd  Vu»  X 

1710  END    IF 

1720  ASSIGN  ?File  TO  D_fileS 

1730  IF  Im=l  AND  Jt=0  THEN 

1740!  READ  DATA  THROUGH  THE  DATA  ACQUISITION  SYSTEM 

1750!  IF  THE  INPUT  MODE  dm)  =  1 

1760  BEEP 

1770  INPUT  "ENTER  FLOWMETER  READING". Fm 

'1780  OUTPUT  709;"AR  AF60  AL63" 

1790  OUTPUT  709:  "AS  SA" 

1800  EtD=0 

1810  FOR  1=1  TO  20 

1820  ENTER  709:Et 

1830  Etp=Etp+Et 

1840  NEXT  I 

1350  Etp=Etp/20 

1860  OUTPUT  709: "AS  SA" 

1870  Ptran=0 

1880  FOR  1=1  TO  50 

1890  ENTER  709:Pt 

1900  Ptran=Ptran+Pt 

1910  NEXT  I 

1920  Ptran=Ptran/50 

1930  OUTPUT  709: "AS  SA"  , 

1940  ENTER  709:Bvol 

1950  OUTPUT  709: "AS  SA" 

1960  ENTER  709:Bamp 

1970  OUTPUT  709:"AR  AF20  AL24" 

1980  FOR  1=0  TO  4 

1990  OUTPUT  709: "AS  SA" 

2000  ENTER  709:Emf(I) 

2010  Emf (I)=ABS(Emf (I)) 

2020  NEXT  I 

2030  OUTPUT  713:"T1R2E" 

2040  WAIT  2 

2050  ENTER   713:Tn 

2060  OUTPUT  713;"T2R2E" 

2070  WAIT  2 

2080  ENTER  713;T2 

2090  OUTPUT  713:"T1R2E" 

2100  WAIT  2 

2110  ENTER  713:T12 

2120  Tl=(T11+T12)».5 

2130  CLEAR  713 

2140  ELSE 

2150!  READ  DATA  FROM  A  USER-SPECIFIED  FILE  IF  INPUT  MODE  (Im)  =2 

2160  ENTER  @Fi le; Bvol .Bamp .Ptran .Etp .Emf (•) .Fm , T1  , T2 

2170  END  IF 

2180  Tsat=FNTvsv((En)f^O)+Emf  (1  ))».5) 
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2190  Ti=FNT./sv(Enif  <2)) 

2200  Grad  =  FNGfad<(Tl-fT2)*.5) 

2210  To=Ti+ABS(Etp)/(lO-Grad)»1 .E+6 

2220  IF   Jj=0    THEN 

2230  Er1=ABS<Ti-T1) 

2240  PRINTER   IS    1 

2250  PRINT  USING   T1  =   '•"  .DD  .  3D"  :  T1 

2260  PRINT   USING   Ti  =    •"•  .DD .DD" : Ti 

2270  IF   Erl>.5    THEN 

2280  BEEP 

2290  PRINT  "OCT  AND  TC  DIFFER  MORE  THAN  0.5  C" 

2300  BEEP 

2310  INPUT  "OK  TO  GO  AHEAD  (  I  =Y  .0  =  N)7"  ,Olcl 

2320  END  IF 

2330  PRINT  USING  DT  (OCT)    =  "" .Z .3D" : T2-T1 

2340  PRINT  USING  DT  (T-PILE)  =  "" .Z .3D" : To-Ti 

2350  IF  Ok  1=0  AND  Erl>.5  THEN  3600 

2360  Er2=ABS((T2-T1)-<To-Tx))/(T2-T1) 

2370  IF  Er2>.05  THEN 

2380  BEEP 

2390  PRINT  "OCT  AND  T-PILE  DIFFER  MORE  THAN  57." 

2400  BEEP 

2410  INPUT  "OK  TO  GO  AHEAD  <  1  =Y .0  =  N)7" .0k2 

2420  IF  0k2=0  AND  Er2>.05  THEN  3600 

2430  END  IF 

2440  PRINTER  IS  701 

2450  END  IF 

2460!  CALCULATE  THE  LOG-MEAN-TEMPERATURE  DIFFERENCE 

2470  IF  Itm=1  THEN 

2480  Tf=T1 

2490  T1=T2 

2500  END  IF 

2510  IF  Itm=2  THEN 

2520  Tf=Ti 

2530  Tl=To 

2540  END  IF  ' 

2550  IF  Itm=3  THEN 

2560  Tf  =  <TH-Ti)».5 

2570  TI=(T2+To)*.5 

2580  END  IF 

2590  Tavg=(Tf+Tl)-.5 

2600  Trise=Tl-Tf 

2610  LfTitd=Trise/LQG<(Tsat-Tf  )/(Tsat-TI)> 

2620  Cpu=FNCpw<Tavg) 

2630  Rhow=FNRhou(Tavg> 

2640  Ku)=FNKu(Tavg) 

2650  Muuia=FNMiiy(  Tavg) 

2660  Prij=FNPrw<Tavg) 

2670  Mdt  =  5.00049E-3+6.9861937E-3»Ftn 

2680  Md=Mdt» <1 .0365-Tf»(1 .96644E-3-Tf»5 .252E-6) )/ .995434 

2590  Vf=Md/Rhou 

2700  Vu=Vf/<PI»Dr2/4) 

2710  Trise  =  Trise-.004*'-/u>'2 

2720  Q=Md»CpwTrise 

2730  Qp=0/(PI»Oo»L) 

2740  Uo=Op/Lmtd 

2750  Re=Rhotn*Vu»Di/Muuia 

2760  Fe1=0 

2770  Fe2=0 

2780  Cf=l 

2790  Tu>o=Tsat-5 
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2300  Tf  ilm  =  Tsat/2  +  T.jo/2 

2810  Kf=FNKui<Tf  ilm) 

2820  Rhof=FNRhou(Tf ilm) 

2830  Muf=FNMau(Tf ilm) 

2840  Hfgp=FNHfg(Tsat>-t-.68»FNCpu<Tf  ilm)»<Tsat-Tijo) 

2850  Netj=Kf*(Rhof  *2»9.799»Hfgp/(Muf*Do*0p))  ^3333 

2860  Ho=.655«Neuj 

2870  Tuoc=Tsat-Qp/Ho 

2880  IF  ABS((Tii»oc-Tuo)/Tu»oc)>.001  THEN 

2890  Tuo=Tuoc 

2900  GOTO  2800 

2910  END  IF 

2920  Cf=1.0 

2930  Q(nega=Re'.S»Prij~.3333*Cf 

2940  Hi  =  Kuj/Di*Ci»Gmeqa 

2950  IF  lfe=0  THEN  3040 

2960  P1=PI»<Di+Do) 

2970  P2=PI*(Di+Dr) 

2980  ft1=<Do-Di)»PI*<Di+Do)*.5 

2990  fi2=(Dr-Di)*PI*(Di+Dr)».5 

3000  M1=(Hi*P1/<Kca»A1 ))'.5 

3010  f12=(Hi*P2/(Kcu*A2))\5 

3020  Fe1=FNTanh(M1-L1)/(M1»L1 > 

3030  Fe2=FNTanh(t12»L2)/(f12-L2) 

3040  Dt=Q/(PI»Di-<L+L1»Fe1+L2»Fe2)*Hi) 

3050  Cfc=(Muua/FNMau(Tavg+Dt))^.14 

3060  IF  ABS(<Cfc-Cf )/Cfc)>.D1  THEN 

3070  Cf=(Cf+Cfc)*.5 

3080  GOTO  2930 

3090  END  IF 

3100  X  =  Do*Neu/<Omega*i<u) 

3110  Y=Neu.*(1/Uo-Rm) 

3120!  COMPUTE  COEFFICIENTS  FOR  THE  LEAST-SOUflRES-FIT  STRAIGHT  LINE 

3130  IF  Jp=1  THEN  OUTPUT  ^Filep:X.Y 

3140  Sx=Sx+X 

3150  Sy=Sy+Y  ' 

3160  Sxs=Sxs+X»X 

3170  Sxy=Sxy+X*Y 

3180?  STORE  RAW  DATA  IN  A  USER-SPECIFIED  FILE  IF  INPUT  MODE  (Im>  =  1 

3190  IF  Im=1  AND  Jj  =  0  THEN  OUTPUT  <«Fi  ie:  Bvol  .Bajnp  .Ptran  .Etp  .Enf  (  »)  ,Fm.  T1  ,  T2 

3200  IF  Jj=0  OR  Jp=1  THEN  PRINT  USING  "aX .5(2X .3D .DD) .2(2X .D .50)" : Tf , TI . Tsat .Lm 

td.Vu.X.Y 

3210  BEEP 

3220  J=J+1 

3230  IF  Im=1  AND  Jj=0  THEN 

3240  INPUT  "DO  YOU  HAVE  MORE  DATA  < 1 =Y . 0=N)7" ,Go_on 

3250  Nran=J 

3260  IF  Go_on=1  THEN  1730 

3270  ELSE 

3280  IF  J<Nrun  THEN  1730     - 

3290  END  IF 

3300  Sl  =  (Nrun*Sxy-SySx)/(Nnin»Sxs-Sx  2) 

3310  Ac=(Sy-Sl»Sx)/Nr.in 

3320  Cic=1/SI 

3330  Jj=Jj+1 

3340  IF  Jp=1  THEN  Jp=2 

3350  IF  ABS(<Cic-Ci)/Cic)>.001  THEN 

3360  Ci=(Cic+Ci)».5 

3370  PRINT  USING  "10X .""Tntermediate  Sieder-Tate  coefft  =  ■•" . Z .40" :Ci 

3330  GOTO  1600 
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3^00  IF  Jp  =  0  THEM  .Jp=l 

3410  END  IF 

3420  IF  Jp=l  THEN  1600 

3430  Ci=(Ci+Cic)*.5 

3440  PRINT 

3450  PRINT  USING  "1  OX . ""Sieder-Tate  coefficient         =  "" . Z . 4D" :Ci 

3460  PRINT 

3470  PRINT  USING  "1  OX. '•"Least-Squares  Line: 

3480  PRINT  USING  "lOX,""   Slope     =   "" . Z .5DE . " :SI 

34S0  PRINT  USING  •MOX.""'   Intercept  =  •'•'.MZ.SDE  .  "  :  Ac 

3500  PRINT 

3510  IF  Ini=l  THEN 

3520  BEEP 

3530  PRINT  USING  "lOX , ""NOTE :  "",ZZ.""  data  runs  are  stored  in  file  "".8fl":J,D 

f  ileS 

3540  ELSE 

3550  PRINT  USING  "10X,""NOTE:  Above  analysis  yas  performed  for  data  in  file  "", 

10A":D_f ileS 

3560  END  IF 

3570  PRINT  USING  "15X.""Plot  data  are  stored  in  file  "" , 1 Ofl" ;Plot5 

3580  ASSIGN  <^File  TO  * 

3590  ASSIGN  SFileo  TO  * 

3600  END 

3610  DEF  FNRhou»<T) 

3620  Ro=I006.35724-T»( . 774489-T*<2 .262459E-2-T*3 . 03304E-4> ) 

3630  RETURN  Ro 

3640  FNEND 

3650  DEF  FNPru>(T) 

3660  Pr'j=FNCDu<T)*FNMuu(T)/FNKu»(T) 

3670  RETURN  Prw 

3680  FNEND 

3690  DEF  FNMuu(T> 

3700  A=247. 3/(1+133. 15) 

3710  Mu=2.4E-5*10"A 

3720  RETURN  Mu  " 

3730  FNEND 

3740  DEF  FNKu(T) 

3750  X=(T+273. 15)7273.15 

3760  Ku>=-.92247+X*(2.3395-X»<  1  .8007-X»<  .52577- .  07344*X) )  > 

37-70  RETURN  Kw 

3780  FNEND 

3790  DEF  FNTvsv(Emf) 

3800  COM  /Cc/  C(7) 

3810  Suin=C<0) 

3820  FOR  1=1  TO  7 

3830  Sum  =  Su;n+C(I)»Emf'I 

3840  NEXT  I 

3850  RETURN  Sum 

3860  FNEND 

3870  DEF  FNCpu<T) 

3380  Cpu=(4.2n20858-T»(2.26826E-3-T»<4.42351E-5+2.71428E-7)))*I000 

3890  RETURN  Cpw 

3900  FNEND 

3910  DEF  FNTanh(X) 

3920  P=EXP(X) 

3930  Q=EXP<-X) 

3940  Tanh=(P-t-0)/(P-Q) 

3950  RETURN  Tanh 

3960  FNEND 

3970  DEF  FNGrad(T) 
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3980  COM  /Cc/  C<7) 

3990  Grad=37.9853+.104388*T 

4000  RETURN  Grad 

4010  FNEND 

4020  DEF  FNFvstd) 

4030  F=466.444+T»<7. 09451-1*1 .S5808E-2) 

4040  RETURN  F 

4050  FNEND 

40S0  DEF  FNHfg(T) 

4070  Hfg=2477200-2450»(T-10) 

4080  RETURN  Hfg 

4090  FNEND 
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PROGRAM    DRP 


1000!  FILE   NPtlE:    DRP 

1010?  REVISED:   November  18.  1933 

1020! 

1030  COM  /Cc/  C(7) 

1040  DIM  Emf(10> 

1050  DATA  0. 10086091 .25727.94369.-767345.3295.78025535.31 

1060  DATA  -9247486589,6.97688E+1 I .-2.66192E+13.3.94078E+14 

1070  READ  C(») 

1080  Di=.0127      •  Inside  diameter  of  test  tube 

1090  Do=. 01905     ?  Outside  diameter  of  test  tube 

1100  Dr=. 015875    •  Outside  diameter  of  the  outlet  end 

1110  Dssp=.1524    ?  Inside  diameter  of  stainless  steel  test  section 

1120  Ax=PI^Dssp^2/4-PI*Do-L 

1130  L=. 130175     !  Condensing  length 

1140  L1  =  . 060325    !  Inlet  end  "fin  length"' 

1150  L2=. 034925    ?  Outlet  end  "fm  length" 

1160  Kcu=385       ?  Thernal  conductivity  of  Copper 

1170  Ci=.034      !  Sieder-Tate  coefficient 

1180  Rm=Do*L0G<Do/Di>/(2*Kcu>  !  Wail  resistance  based  on  outside  area 

1190  PRINTER  IS  701 

1200  CLEAR  709 

1210  BEEP 

1220  INPUT  "ENTER  MONTH.  DATE  AND  TIME  (MM:DD:HH:MM:SS)" ,Date5 

1230  OUTPUT  709:"TD";DateS 

1240  OUTPUT  709;"TD" 

1250  ENTER  709:DateS 

1260  PRINT  "  Month,  date  and  time  :";Date$ 

1270  PRINT 

1230  PRINT  USING  "1  OX. ""NOTE:  Program  name  :  DRP 

1290  BEEP 

1300  INPUT  "ENTER  DISK  NUMBER". Dn 

1310  PRINT  USING  "15X.""Disk  numcer   =  "",DD":Dn 

1320  BEEP 

1330  INPUT  "ENTER  INPUT  MODE  ( 1 =3054A,2=FILE)". Im 

1340  IF  Im=1  THEN 

1350  BEEP 

1360  INPUT  "GIVE  A  NAME  FOR  THE  RAW  DATA  FILE" ,D_f ile$ 

1370  CREATE  BDAT  D_fileS.15 

1380  ASSIGN  5File  TO  D_fileS 

1390  BEEP 

1400  INPUT  "ENTER  GEOMETRY  CODE  < 1 =FINNED  .  0=PLAIN)"  ,  Ifg 

1410  OUTPUT  §File:Ifg 

1420  IF  lfg=0  THEN 

1430  BEEP 

1440  INPUT  "WALL  TEMPERATURE  MEASUREMENT  ( 1 =Y ,0=N>" . lut 

1450  ELSE 

M60  BEEP 

1470  INPUT  "ENTER  FIN  PITCH,  WIDTH  AND  HEIGHT" .Fp.Fu .Fh 

1480  END  IF 

M90  IF  lfg=0  THEN  OUTPUT  SFiletlut 

1500  IF  Ifg=1  THEN  OUTPUT  SFile :Fp,Fu .Fh 

1510  ELSE 

1520  BEEP 

1530  INPUT  "GIVE  THE  NAME  OF  THE  EXISTING  DATA  FILE" .D_f i ieS 

1540  PRINT  USING  "15X.""This  analysis  was  performed  for  data  in  file  "".10A";D 
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1550  BEEP 

1560  INPUT  "ENTER  THE  NUMBER  OF  RUNS  STORED" .Nrun 

1570  ASSIGN  'aPiie  TO  D_file£ 

1580  ENTER  3File:Ifg 

1590  IF  lfg=0  THEN  ENTER  «File:Iijt 

1B00  IF  Ifg=l  THEN  ENTER  iFi  ie:Fp  .F.j.Fh 

1610  END  IF 

1520  IF  lfg=0  THEN 

1530  BEEP 

1640  INPUT  "WANT  TO  CREATE  A  FILE  FOR  Nr  vs  F  ( 1 =Y . 0=N)7" . Inf 

1650  ELSE 

1660  lnf=0 

1670  END  IF 

1680  IF  Inf=1  THEN 

1690  BEEP 

1700  INPUT  "GIVE  A  NAME  FOR  Nr  vs  F  FILE".NrfS 

1710  CREATE  BOAT  NrfS.2 

1720  ASSIGN  BNrf  TO  NrfS 

1730  END  IF 

1740  BEEP 

1750  INPUT  "ENTER  OPTION  ( 1 =0CT.2=T-PILE.3=AWE)" . Itn 

1760  BEEP 

1770  INPUT  "ENTER  OPTION  FOR  END-FIN  EFFECT  (I =Y .0=N)"  .  Ife 

1780  IF  Itm=1  THEN  PRINT  USING  "16X.""This  analysis  uses  OCT  readings 

1790  IF  Itm=2  THEN  PRINT  USING  "16X.""This  analysis  uses  T-PILE  readings 

1800  IF  Itm=3  THEN  PRINT  USING  "16X.""This  analysis  uses  average  of  QCT  and  T-P 

ILE  readings 

1810  IF  Ife=1  THEN  PRINT  USING  "16X.""Thi3  analysis  inclLides  end-fin  effect 

1820  IF  Ife=n  THEN  PRINT  USING  "16X.""This  analysis  neglects  end-fin  effect 

1830  PRINT  USING  "16X .""Sieder-Tate  coefficient   =  "".Z.4D":Ci 

1840  BEEP 

1850  INPUT  "GIVE  A  NAME  FOR  PLOT  DATA  FILE" .P_f ileS 

1860  CREATE  BOAT  P_fileS.5 

1870  ASSIGN  *Filep  TO  P_fileS 

1880  IF  Iut=1  THEN 

1890  BEEP 

1900  INPUT  "GIVE  A  NAME  FOR  WALL  TEMPERATURE  FILE",;4tfS 

1910  CREATE  BOAT  WtfS.5 

1920  ASSIGN  itFilel  TO  WtfS 

1930  END  IF 

1940  BEEP 

1950  INPUT  "ENTER  OUTPUT  VERSION  < T=SHORT .2=L0NG)" . lov 

1960?  IF  Im=1  THEN 

1970?  OUTPUT  iFile;Ifg 

1980?  IF  lfg  =  0  THEN  OUTPUT  «File:Iujt 

1990?  IF  Ifg=I  THEN  OUTPUT  'aFi  le:Fp  .Fuj.Fh 

2000?  ELSE 

2010?  ENTER  ?Fiie:Ifg 

2020?  IF  lfg  =  0  THEN  ENTER  («File;Iut 

2030?  IF  Ifg=I  THEN  ENTER  ^^FileiFp  .Fui  .Fh 

2040?  END  IF 

2050  IF  lfg=0  THEN 

2060  PRINT  USING  "16X,""Tube  type     :  PLAIN 

2070  ELSE 

2080  PRINT  USING  "16X.""Tuhe  type     :  FINNED 

2090  PRINT  USING  "16X.""Fin  pitch,  uiidth.  and  height  (mm):  ""  .DD  .  D  ,2X  .  Z  .DO  .2X  .  Z 

.  DD"  ;Fp  ,Fuj  .Fh 

2100  END  IF 

2110  J=0 

2120  IF  Iov=I  THEN 
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2130      PRINT 

2!  40      IF    Inf=I    THE^J 

2!  50  PRINT  USING  "1  OX.  ""Data  ^>u      Uo       Ho        Qp      Vv    F     Hr 


2''SQ  PRINT  USING  "tOX.""  #   (n/s)  (W/n^2-<)(W/r»  2-0   (W/m^2)  <m/s) 

2170  ELSE 

2t30  PRINT  USING  "TOX. ""Data  Vu        Uo  Ho  Qp       Vv 

2130  PRINT  USING  "10X.""  #   (m/s)    (W/m"2-K)     <W/m  2-K)     (W/m  2)    (tn/s)"" 

2200  END  IF 

2210  END  IF 

2220  Sx=0 

2230  Sy=0 

2240  Sxs=0 

2250  Sxy=0 

2260  Repeat;! 

2270  J=J+1 

2280  IF  Im=I  THEN 

2230  BEEP 

.2300  INPUT  "LIKE  TO  CHECK  NG  CONCENTRATION  (I =Y . 0=N)?" .Ng 

2310  BEEP 

2320  INPUT  "ENTER  FLOWMETER  READING". Fm 

•2330  OUTPUT  703;"AR  AFSO  ALS3  VR5" 

2340  OUTPUT  709: "AS  SA" 

2350  ENTER  709:Etp 

2360  OUTPUT  709: "AS  SA" 

2370  Vtran=0 

2380  FOR  1=1  TO  50 

2390  ENTER  709 :Vt 

2400  Vtran=Vtran+Vt 

2410  NEXT  I 

2420  Vtran=Vtran/50 

2430  OUTPUT  709; "AS  SA" 

2440  ENTER  709:Bvol 

2450  OUTPUT  709: "AS  SA" 

2460  ENTER  709;Bamp 

2470  IF  lut=0  THEN  OUTPUT  709;"AR  AF20  AL24  VR1" 

2480  IF  Iwt=1  THEN  OUTPUT  709;"AR  AF20  AL30  VR1" 

2490  IF  luit  =  0  THEN  Nn=4 

2500  IF  Iut  =  1  THEN  Nti  =  10 

2510  FOR  1=0  TO  Nn 

2520  OUTPUT  709; "AS  SA" 

2530  IF  I>4  THEN 

2540  Se=0 

2550  FOR  K=0  TO  10 

2560  ENTER  709 ;E 

2570  Se=Se+E 

2580  NEXT  K 

2590  Etnf(I)=ABS(Se/10) 

2600  ELSE 

2610  ENTER  709:E 

2620  Emf (I)=ABS<E) 

2630  END  IF 

2640  NEXT  I 

2650  OUTPUT  713;"TtR2E" 

2660  WAIT  2 

2570  ENTER  713:Tn 

2680  OUTPUT  713;"T2R2E" 

2590  WAIT  2 

2700  ENTER  713;T2 
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^/lU  uuiPui  /Io;"riK2E" 

2720  WAIT  2 

2730  ENTER  713:T12 

27A0  T1=<T1 1+T12>».5 

2750  IF  Ng=0  THEN  2800 

27S0  BEEP 

2770  INPUT  "ENTER  MANOMETER  READING  (HL  .HR  .HRH)"  .HI  ,Hr  .Hr-j 

2780  Phg=Hi+Hr 

2790  Puater^Hr-Hrui 

2800  ELSE 

2810  IF  Ifg=1  OR  lut=0  THEN 

2820  ENTER  ^File :Bvol .BamP.^tran.Etp .Emf (0) .Emf ( ! ) .Emf (2) ,Emf <3) .Emf ( 4 ) .Fm . Tt . T 

2.Phg.Puater 

2830  END  IF 

2840  IF    lfg=0   AND    Iijt=l    THEN   ENTER   '?>Fi  ie  :  Bvoi  .Bamp  .  Vtran.Etp  .Emf  («  >  .Fn  .  TI  .  T2  .Ph 

2850  IF   J=1    OR   J=10   OR   J=20    OR   J=Nrun    THEN 

2860  Ng=1 

2870  ELSE 

•2880  Ng=0 

2890  END    IF 

2900  END    IF 

2910  TsteajT>  =  FNTvsv((En»f(0)+Enf{1))*.5)       !    COMPUTE   STEAM    TEMPERATURE 

2920  Troom=FNTvsv<Emf (3)) 

2930  IF    Iu.t=1    THEN 

2940  Ttjm=0. 

2950  FOR   1=0   TO   5 

2960  Tu(I)=FNTvsv<Efnf  (1+5)) 

2970  Tum=Tu»m+Tu(I) 

2980  NEXT    I 

2990  Tujm=Tu»m/6 

3000  END    IF 

3010  Tcon=FNTvsv(Emf (4)) 

3020  Psat=FNPvst(Tsteam> 

3030  Rohg=I3529-I22»<Troom-26.85)/50 

3040  Rouater=FNRhou(Troom) 

3050  Ptest  =  (Phg*Rohg-Pii>ater*Rouater)*9.  799/ 1000 

3060  Pmm=Ptest/133.322 

3070  Plcm=Ptesf  1  .E-3 

3080  Plcs=Psaf  1  .E-3 

3090  Pkt=FNPvsv(Vtran)»1 .E-3 

3100  Tsat=FNTvsp(Ptest) 

3110  Vst=FNVvst(Tsteajn) 

3120  Ppng=(Ptest-Psat)/Ptest 

3130  Ppst=l-Ppng 

3140  Mfng=l/(1+18. 015/28. 97-Psat/<Ptest-Psat)) 

3150  Vfng=Mfng/(1 .608-.D08»Mfng) 

3160  Mfng=Mfng-!00 

3170  Vfng=Vfng*100 

3180  BEEP 

3190  IF  Iov=2  THEN 

3200  PRINT 

3210  PRINT  USING  "1  OX . ""Data  set  number  =  "".DD";J 

3220  PRINT 

3230  END  IF 

3240  IF  Iov=2  AND  Ng=l  THEN 

3250  PRINT  USING  "10X.""       P         Psat    Ptran   Tmeas    Tsat        N 

G  •/. 

3260  PRINT  USING  "10X.""  (mm)     (kPa>    (kPa)    (kPa.)  <C)     (C)    Molal 

Mass 

3270  PRINT  USING  "1  OX .5(3D .DD . 2X) .2(3D .DD .2X) .2(M3D . D .2X>" ;Pmn .okm .Pks .Pk t . Tste 
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an ,  Ts 

at . Vf ng .Mf ng 

3280 

PRINT 

3290 

END  IF 

3300 

IF  f1fng>.5  THEN 

3310 

BEEP 

3320 

PRINT 

3330 

IF  Itn=1  THEN 

3340 

BEEP 

3350 

PRINT 

3360 

PRINT  USING  "1  OX. ""Energize  the  vacuum  system  

3370 

BEEP 

3380 

INPUT  "QK  TO  ACCEPT  THIS  RUN  ( 1  =Y  ,0=N)'?"  ,0k 

3390 

IF  Ot=0  THEN 

3400 

BEEP 

3410 

DISP  "NOTE;  THIS  DATA  SET  WILL  BE  DISCARDED??  " 

3420 

WAIT  5 

3430 

GOTO  2280 

3440 

END  IF 

3450 

END  IF 

3460 

END  IF 

3470 

IF  Im=I  THEN 

3480 

IF  Ifg=1  OR  lut=0  THEN 

3490 

OUTPUT  ■^File;Bvol,Bamp,Vtran.Etp.Emf<0).Emfn>.Emf(2) 

.Emf<3), 

Emf(4) .Fn 

..TI, 

T2.Phg.Puater 

3500 

END  IF 

3510 

IF  lfg=0  AND  Iu>t  =  1  THEN  OUTPUT  ^F  i  Ie:Bvoi  .Bamp  .^^tran  . 

Etp,Emf( 

•>.Frr>,T1  . 

T2.P 

hg.Puater 

3520 

END  IF 

3530 

IF  lfg=0  AND  Iuit=1  THEN  OUTPUT  iFi  le' ;  Tu»( ») 

3540! 

ANALYSIS  BEGINS 

3550 

Ti=FNTvsv(Emf (2)) 

3560 

Grad=FNGrad((T1+T2)-.5) 

3570 

To=Ti+ABS<Etp)/<10»Grad)*1 .E+6 

3580 

Er1=ABS<Ti-T1) 

3590 

PRINTER  IS  1 

3600 

PRINT  USING  T1  (OCT)    =  "".DD.aD":T1 

3610 

PRINT  USING  Ti  (TO     =  "".DD.3D":Ti 

3620 

IF  Er5>.5  THEN 

3630 

BEEP 

3640 

PRINT  "QCT  AND  TC  DIFFER  BY  MORE  THAN  0.5  C" 

3650 

BEEP 

3660 

INPUT  "QK  TO  GO  AHEAD  < 1 =Y .a=N)7" ,0k  I 

3670 

END  IF 

3680 

PRINT  USING  DT  <GCT)    =  "" . Z .3D" : T2-T1 

3590 

PRINT  USING  DT  (T-PILE)  =  "" . Z .3D" : To-Ti 

3700 

IF  Gkl=0  AND  Er1>.5  THEN  5100 

3710 

Er2=ABS<(T2-T1 )-( To-T i ) )/( T2-T1  ) 

3720 

IF  Er2>.05  THEN 

3730 

BEEP 

3740 

PRINT  "QCT  AND  T-PILE  DIFFER  BY  MORE  THAN  57." 

3750 

BEEP 

3760 

INPUT  "OK  TO  GO  AHEAD  <  1  =Y  .0=N)'7"  ,0k2 

3770 

IF  0k2=0  AND  Er2>.05  THEN  5100 

3780 

END  IF 

3790 

PRINTER  IS  701 

3800 

IF  Itm=1  THEN 

3810 

T1i=T1 

3820 

T2o=T2 

3830 

END  IF 

3840 

IF  Itm=2  THEN 

3850 

T1i=Ti 
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3860  T2o=To 

3870  END  IF 

3880  IF  Itm=3  THEM 

3890  T1 i=(T1+Ti>*.5 

3900  T2o=(T2+To)*.5 

3910  END  IF 

3920  Tavg=(T:i+T2o)».5 

3930  Cpij=f^NCDy(Tavg) 

3340  Rhou=FNRhou<Tavg) 

3950  Md=5.00049E-3+6.9861937E-3»Fm 

3960  f1d=Md»( 1 .0365-1 .96b44E-3»Tavg+5 .252E-6-Tavg  2)/ .995434 

3970  Mf=Md/Rhou 

3980  y.j=Hf/(PI»Di'2/4) 

3990  T2o=T2o-( .0138+.001»^u*2) 

4000  Q=Md*Cpuc<T2o-Tl i) 

4010  Qp=Q/(PI*Do»L) 

4020  Ku  =  FNKui(Tavg) 

4030  MiiLj=FNMuu(Tavg) 

4040  Rei=Rhoii)*Vij*Oi/riuuj 

4050  Prui=FNPruj(Tavg) 

4060  Fe1=0. 

4070  Fe2=0. 

4080  Cf=!. 

4090  Hi=Ku.*Ci/Di»Rei^.8»Pru\3333»Cf 

4100  Dt=Q/(PI*Di»<L+L1*Fe1+L2»Fe2)*Hi) 

41  10  Cfc=<Maij/FNMuu(Tavg+Dt) )  M4 

4120  IF  ABS((Cfc-Cf)/Cfc)>.01  THEN 

4130  Cf=(Cf+Cfc)».5 

4140  GOTO  4090 

4150  END  IF 

4150  IF  lfe=0  THEN  GOTO  4250 

4170  P1=PI*(Di+Do) 

4180  A1=<Do-Di)»PI»^Di+Do)».5 

4190  M1=<Hi»P1/(Kcu»A1 )>* .5 

4200  P2=PI*<Di+Dr) 

4210  fl2=<Dr-Di)*PI-(DL+Dr)*.5 

4L20  M2=(Hi»P2/(Kca*A2>)  .5 

4230  F8l=FNTanh(M1-L1 )/(M1»L1 ) 

4240  Fe2=FNTanh(ri2»L2)/(M2*L2) 

4250  Lmtd=<T2o-T1  i  )/LOG<  <  Tsteam-TI  i  )/(Tsteaji>-T2o)  > 

4260  Uo=Q/(Lmtd»PI-Oo*L) 

4270  Ho= 1 / ( 1 /Uo-Do*L/ ( D i-( L>L I -Fe 1 +L2«Fe2 ) »H i ) -Rm) 

4280  Dtc=Q/(PI»Di»(L+Ll»Fel+L2»Fe2>»Hi) 

4290  IF  ABS((Dtc-Dt)/Dtc)>.01  THEN  4090 

4300  Hfg=FNHfg<Tsteam) 

4310  Tuo  =  TsteajTi-Qp/Ho 

4320  Tfiim=Tsteajn/3+2*Tyo/3 

4330  Kf=FNKb>(Tf  iltn) 

4340  Rhof=FNRhou(Tf ilm) 

4350  Muf=FNMuu(Tfiini) 

4360  Hpq=.651»Kf«<Rhof*2*9.81*Hfg/(Muf*Do*Qp)) ".3333 

4370  Y=Hpg»Qp-.3333 

4380  X=Qp 

4390  Sx=Sx+X 

4400  Sy=Sy+Y 

4410  Sxs=Sxs+X'2 

4420  Sxy=Sxy+X»Y 

4430  OUTPUT  §Filep:Qp.Ho 

4440  01=500   ?  TO  BE  MODIFIED 

4450  Qloss=Ql/(100-25)»(Tsteajn-TrooiD)   ?  TO  BE  MODIFIED 

4460  Hfc=FNHf(Tcon) 
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4470  Hf=FNHf (Tstean) 

4^80  Mdv=0 

4490!  Bp=(Bvo1*100) ^2/5.7S   !  BOILER  POUER  IM  Watts 

4500  Bp=(8voi*100)'2/5.7S 

4510  Mdvc=(  (Bp-Qloss)-t1dv*(Hf-rifc)  )/Hfa 

4520  IF  ABS(<Mdv-Mdvc)/Mdvc)>.01  THEN 

4530  f1dv=<lidv+Hdvc>*.5 

4540  GOTO  4510 

4550  END  IF 

4560  Mdv=(Mdv+Mdvc)».5 

4570  Vg=FNVvst(Tsteam) 

4580  Vv=Mdv*Vg/Ax 

4590  IF  Inf=I  THEN 

4G00  F=(9.799»Do*Muf*Hfg)/(Vv'2»Kf*<Tsteam-Tajo)) 

4610  Nu=Ho*Do/Kf 

4620  Ret=Vv*Rhof*Do/Muf 

4630  Nr=Na/Ret^.5 

4640  END  IF 

4650  IF  Inf=I  THEN  OUTPUT  aNrf;F.Nr 

4650  IF  Iov=2  THEN 

4670  PRINT  USING  •'lOX.""   T  (Inlet)      Delta-T 

4680  PRINT  USING  "10X.""  QCT     TC    OCT  T-PILE 

4690  PRINT  USING  "1  OX .2(DD .DD .2X) .2(Z .3D .2X)" ; T1 . Ti . T2-T1 . To-T i 

4700  PRINT  USING  "10X.""  Vu      Rei         Hi         Uo         Ho         q 

Vv 

47!0  PRINT  USING  "lOX.  Z  .DD  .  1  X  .5<MZ  .  3DE  .  IX)  .MZ  .DD"  ;U(j.Rei  .Hi  .Uo  .Ho  .Qp  .  Vv 

4720  END  IF 

4730  IF  Iov=I  THEN 

4740  IF  Inf=1  THEN 

4750  PRINT  USING  "1 IX  .DD  .2X  ,Z  .DD  .2X  ,2<5D  .  D  .2X)  .  Z  .3DE  .  IX  .  Z  .DD  .2<  1  X  .3D  .DD)"  :  J  .  ^;.j, 

Uo.Ho.Op.Vv.F.Nr 

4760  ELSE 

4770  PRINT  USING  "MIX  .DD  .2X  .Z  .DD  .2X  .2(f1D  .  4DE  .2X)  .  Z  .  3DE  .3X  .  Z  .  DD"  :  J  .  Vu.Uo  .Ho.Op  .  V 

4780  END  IF 

4790  END  IF  ^ 

4800  IF  Im=1  THEN 

4810  BEEP 

4820  INPUT  "WILL  THERE  BE  ANOTHER  RUN  H =Y .0=N)7" .Go_on 

4830  Nrun=J 

4840  IF  Go_on=I  THEN  Repeat 

4850  ELSE 

4860  IF  J<Nriin  THEN  Repeat 

4870  END  IF 

4880  IF  lfg=0  THEN 

4890  PRINT 

4900  Sl  =  (NruTTSxy-Sy*Sx)/<NruiTSxs-Sx'2) 

4910  Ac=<Sy-Sl*Sx)/Nran 

4920  PRINT  USING  "t OX. ""Least-Squares  Line  for  Hna  vs  q  curve: 

4930  PRINT  USING  "10X.""   Slope      =  "" .MD .4DE" ;SI 

4940  PRINT  USING  "10X,""   Intercept   =  "" ,MD . 4DE" : Ac 

4950  END  IF 

4960  IF  Im=1  THEN 

4970  BEEP 

4980  PRINT 

4990  PRINT  USING  "1  OX .""NOTE :  "".II.""   data  runs  yere  stored  in  file  "".10A":J. 

D_f ileS 

5000  END  IF 

5010  BEEP 

5020  PRINT 

5030  PRINT  USING  "I0X.""NOTE:  "".ZZ.""  X-Y  pairs  .jere  stored  in  plot  data  file 
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•••'.lOfl 

5040 

5050 

5060 

5070 

5030 

5090 

5100 

5110 

5120 

5130 

5U0 

5150 

51S0 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

5360 

5370 

5380 

5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

5470 

5480 

5490 

5500 

55  !0 

5520 

5530 

5540 

5550! 

5560 

5570 

5580 

5590 

5600! 

5610 

5620 

5630 


pairs 


s  of  Nr-F  are  stored  m  file 


MA":J.:JrfS 


J.P_f lies 
IF  Inf=I  THEN 
PRINT  USING  ••16X.ZZ 
END  IF 

ASSIGN  ?File  TO  • 
ASSIGN  *File1  TO  * 
ASSIGN  ^Fiiep  TO  - 
END 

DEF  FNPvst(Tsteam) 
DIM  K(8) 

DATA  -7 . 691 234564 . -26 . 08023696 . -  1 58 . 1 706546 .64 . 23235504 .-118. 9646225 
DATA  4.16711732.20.9750676,1 .E9.5 
READ  K(-) 

T=<Tsteajii+273. 15)7647.3 
SLun=0 

FOR   N=0    TO   4 

SLUii=Siun+K(N)*(1-T)"<N+n 
NEXT   N 

Br=Sum/<T*<1+K(5>»(1-T)+K(6)*(1-T) "2) )-< 1 -T) /(K< 7)*( 1 -T) -2+K<3) ) 
Pr=EXP(Br) 
P=22120000»Pr 
RETURN  P 
FNEND 

DEF  FNHfg(T) 
Hfg=2477200-2450»<T-I0) 
RETURN  Hfg 
FNEND 

DEF  FNMauj(T) 
A=247.3/(T+133.15) 
Mu=2.4E-5*10"A 
RETURN  Ma 
FNEND 

DEF  FNVvstdt) 
P=FNPvst(Tt) 
T=Tt+273.15 

X=1500/T  ' 

F1=1/(I+T*1 .E-4) 
F2=( 1 -EXP<-X) )*2 .5»EXP(X)/X~ .5 
B=.0015«F1-.000942*F2-.0004882*X 
K=2»P/(461 .52»T) 
V  =  (1+(1+2*B*K)~.5)/'!< 
RETURN  V 
FNEND 

DEF  FNCpuCT) 

Cpu=4.21120858-T»<2.2S826E-3-T*<4.4236IE-5+2.71428E-7«T)) 
RETURN  Cpu-1000 
FNEND 

DEF  FNRhou)(T> 

Ro=999.52946+T»< . 01269-T*(5 .482513E-3-T-1 .234147E-5)) 
RETURN  Ro 
FNEND 

DEF  FNPru)(T) 

Pry=iO  •<  1  .09976605-T»<1  .37'a9326E-2-T»(3  .968875E-5-3.  45026E-7*T) ) ) 
Pru=FNCpw(T)*FNMiiiu(T)/FNK(ii<T) 
RETURN  Pru) 
FNEND 

DEF  FNKu)<T) 

Ku=.5625894+T*<2.2964546E-3-T*(l .509766E-5-4 .0581652E-8»T) ) 
X=(T+273. 15)7273.15 

Ku=-.92247+X»(2.8395-X*( I .8007-X*< . 52577- . 07344»X) ) ) 
RETURN  Kw 
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5640  -MEND 

5G50  DEF  FNTanh(X) 

5650  P=EXP<X) 

5670  0=EXP(-X) 

5S30  Tanh=(P+Q)/(P-Q) 

5690  RETURN  Tanh 

5700  FNEND 

5710  DEF  FNTvsv(V) 

5720  COM  /Cc/  C(7) 

5730  SuLm=C(0) 

5740  FOR  r=!  TO  7 

5750  Sam=Sujn+C(I)»V'r 

5760  NEXT  I 

5770!  T  =  ^;»<  .02635206355-^;»(9.735ni3E-7-V*6.576805E-n)) 

5780  RETURN  Sam 

5790  FNEND 

5300  DEF  FNHf(T) 

5810  Hf=T*(4.203849-T«(5.88132E-4-T*4.55160317E-6)) 

5820  RETURN  Hf*l000 

S830  FNEND 

5840  DEF  FNGrad(T) 

5850  Grad=37.9853+.]04388-T 

5860  RETURN  Grad 

5870  FNEND 

5880  DEF  FNTv3p(P) 

5890  Ta=n0 

5900  Tl=10 

5910  Ta=(Tu-t-Ti)*.5 

5920  Pc=FNPvst(Ta) 

5930  IF  ABS(<P-Pc)/P)>.001  THEN 

5940  IF  Pc<P  THEN  Tl=Ta 

5950  IF  Pc>P  THEN  Ta=Ta 

5960  GOTO  5910 

5970  END  IF 

5980  RETURN  Ta 

5990  FNEND  ' 

6000  DEF  FNPvsv(V) 

6010  P=8133. 5133+2. 23605'E+4*<; 

6020  RETURN  P 

6030  FNEND 
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PROGRAM  TCAL 


'00  !  "ILE  NAME:  TCAL 

110  !  REVISED:    Decemoer  1i.  1983 

!20  ! 

130  COM  /Cr/  C(7) 

UO  DIM  Emf  <10)  .1(10)  .D<10) 

150  DATA  0.10086091.25727.94369.-7573^5.3235.73025595.81 

160  DATA  -9247486589. 6. 97G88En .-2. 55192E13. 3. 94078E1^ 

170  READ  CO 

180  CLEAR  709 

190  BEEP 

200  INPUT  "ENTER  MONTH,  DATE  AND  TIME  (MM:DD  :HH:Mht:SS)"' .BS 

210  J=0 

220  OUTPUT  709:"TD":B3 

230  OUTPUT  709:"TD" 

240  ENTER  709: AS 

250  PRINT   USING   "1  OX  .  ""MoTith  .    date   and    time      =    ■•■M4A";AS 

250  BEEP 

270  INPUT  "ENTER  INPUT  MODE  (I=3054A.  2=FILE)".ini 

280  IF  Im=1  THEN 

290  BEE? 

300  INPUT  "GIVE  A  NAME  FOR  DATA  FILE" .D_f i leS 

310  CREATE  BOAT  D_fiie$.5 

320  ELSE 

330  BEEP 

340  INPUT  "GIVE  NAME  OF  EXISTING  FILE" .D_f i leS 

ocQ  BEEP 

360  INPUT  "ENTER  NUMBER  OF  DATA  RUNS  STORED", Nran 

370  END  IF 

380  BEE? 

390  INPUT  "GIVE  A  NAME  FOR  PLOT  FILE"  .f*_f  i  ieS 

400  CREATE  BDAT  P_fiiaS.5 

410  ASSIGN  iPlot  TO  P_fiieS 

420  ASSIGN  '?iFiie  TO  D_fiieS 

430  IF  I,'T»=1  THEN 

440  BEEP 

450  INPUT  "ENTER  BATH  TEMPERATURE" . T_bath 

460  OUTPUT  709:"AR  AF20  AL30" 

470  FOR  1=0  TO  10 

480  OUTPUT  709: "AS  SA" 

490  ENTER  709:Emf(I) 

500  NEXT  I 

510  OUTPUT  713:"T1R2E" 

520  WAIT  2 

530  ENTER  713:T1 

540  OUTPUT  713:"T2R2E" 

550  WAIT  2 

550  ENTER  713:T2 

570  OUTPUT  ?File:T_bath.Emf(*),T1 .T2 

580  EL  '^F 

590  ENTER  aPi ie: T_bath.Emf ( ») . T1 , T2 

600  END  IF 

b;0  J=J+! 

620  Duja=0 

530  FOR  1=0  TO  10 

640  T(I)=FNTvsv<ABS(E!nt<I))) 

650  D(I)=T_bath-T(I) 

660  IF  I>4  THEN  Duja=Dwa+D<  I) 

670  NEXT  I 

630  Dua=Duja/5 
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690  Dsa=(D<0)+0<  1 ) )*.5 

700  OUTPUT   >5Plot:T_bath.Dsa.Dua 

7!n  PRINT 

720  PRINT  USING  "]OX.""Data  set  number  =  "".DD'^J 

730  PRINT  USING  "1  OX  .  ""Bath  T  (O    QCT-1  (C)    QCT-2  (O 

740  PRINT  USING  "1  OX .3< 3D . 3D . 7X)" ; T_bath . T1 . 12 

750  PRINT  USING  "1  OX  .  ""Thermocouple  readinas  (Deq  O: 

750  PRINT  USING  "t OX .b< 3D .DD . 3X) , 16X" ; T( *) 

770  PRINT  USING  "1  OX . ""Discrepanc les  (Deg  C): 

730  PRINT  USING  "HX  .6<t1Z  .DD  .4X)  ,  i  5X"  :DC) 

790  BEEP 

800  IF  Im=l  THEN 

810  INPUT  "ARE  YOU  TAKING  i^ORE  DATA  ( t  =  r  .  0=N)7"  .Go_on 

820  IF  Go_on=1  THEN  430 

340  IF  J<Nrun  THEN  430 

850  END  IF 

860  PRINT 

870  IF  Im=l  THEM 

380  PRINT  USING  "10X.""NOTE:  "".DD.""  data  sets  are  stored  in  file  "".14A":J.D 

_f  ileS 

390  ELSE 

900  PRINT  USING  "10X.""NOTE:  Above  analysis  ijas  oer  +  ormed  froni  file  "".!4A":D_ 

f  ileS 

9rO  END  IF 

S20  PRINT  USING  "•6X.""Plot  data  are  stored  in  file  "" . 1 OA" : P_f i le5 

930  ASSIGN  'i^ile  TO  * 

940  ASSIGN  !»Plot  TO  •- 

950  END 

960  DEF  FNTvsv<Emf) 

970  COM  /Cc/  C(7) 

980  Suiii=C(0) 

990  FOR  1=1  TO  7 

1000  Sujn=S.jjn+C(I)*Enf  I 

1010  NEXT  I 

1020  RETURN  Sum 

1 030  FNEND 
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